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1.0 INTRODUCTION 


Pulsed plasma thrusters (PFT) using solid teflon propellant have a flight 

demonstrated sijnplicity and reliability and are of increasing interest for 

future flight applications. Earlier versions of these thrusters had total 

impulses of less than 3000 lb -sec and impulse bits of less than 100 plb-sec, 

and hence, were limited to applications such as east -west stationkeeping and 

C21 

fine attitude control of small spacecraft.^ A one-millipound average thrust 
PPT is currently under development to extend tJie capabilities of the pulsed 
plasma thruster to applications on larger spacecraft with longer mission 
duration. ^ This thruster has an impulse bit roughly 50 times larger, and a 
total impulse rouglily 10 times larger, than the earlier versions. In addition, 
its specific impulse, propellant flow rate, efficiency and power are also 
significantly higher. Envisioned applications for this millipound thruster 
include nort}i-sout]i stationkeeping, satellite orbit acquisition and maneu- 
vering, and large structure attitude control. 

Tlie flight experience of the smaller pulsed plasma thrusters has sho\\n 
that the exhaust plume of these thrusters has a negligible effect on spacecraft 
surfaces . The exhaust plume of the one-millipound thruster is of potentially 
greater concern, primarily due to its higher energy and mass. In addition, 
longer duration missions with ever more sensitive instrumentation will aggravate 
my plume contamination problem that may exist. Previous studies have been 
conducted to assess the effect of the one-millipound PPT plume on spacecraft 
surfaces by directly measuring the plume flux towards a spacecraft upstream of 
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the tliruster exhaust plane. Unfortunately, accurate results have been masked 
by a backscattered flux of particles reflected and eroded from the test facility 
vacuum chamber walls. In order to minimize this effect and to develop a more 
accurate measure of the plume -spacecraft interaction, a study was carried out 
at the Jet Propulsion Laboratory using the Molecular Sink Vacuum Facility 
(MOLSINK) , This facility has a gaseous cooled anechoic-type liner 

(MOLTRAP) especially designed to minimize any plume -wall backscatter, thus 
providing aji environment in which accurate plume -spacecraft interaction 
measurements may be made. 

The JPL plume study has been divided into two phases; 

Phase I ; An evaluation of the PPT plume -\i?all backscatter characteristics 
of the MOLSINK facility, a conceptual design for a PPT backflow measurement 
tecfmique, and the development of a low temperature quartz crystal 
microbalance CQCM) desi^ to be used in measuring this backflow. 

Phase II : A study of the plume -spacecraft interaction utilizing the MOLSINK 
facility, including a direct measurement of the plume backflow mass flux 
into the thruster nozzle exit plane at various distances frcmi the thruster 
axis and a measure of the PPT primary plume mass flux profile downstream 
of the thruster nozzle. 

The experiments and analyses of Phase I have been completed and are detailed 
in a previous report. Included is a description of the MOLSINK' facility as 
modified for use with the pulsed plasma thruster, and of the Solar Electric 
Propulsion (SEP) vacuum facility, used for the PPT primary plume studies and 
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preliminary QOl testing. The QCM circuit design and calibration at liquid 
nitrogen tenperatures is also discussed. Backscatter from the MOLTRAP anedioic 
surface was measured in total, and at two specific locations, using QQ! test 
arrays designed for this purpose. The results indicate that the plume wall 
backscatter is highest at the wall areas closest to the thruster axis and 
Cnris to iicgLigiblc values at the plume l)oundary-wall intersection, 4(’ ' oil this 
axis. Tlie total backscatter from the entire MOLTRAP wall area was found to be 
almost 5% of the PPT plune mass. Based on these wall backscatter measurerents, 
a conceptual method was developed for measuring the PPT plume backflow. Because 
of the relatively large plume-wall backscatter in the MOLSINK over the wall 
area directly in the PPT primary plume, an indirect method of measuring the 
backflow is required which avoids measuring the wall backscatter. Figure 1 
shows the conceptual technique, which uses a collimated QCM to make the plume 
backflov; measurement. This collimated QCM is rotated around a fixed point at 
the entrance to the collimator. The collimator aperture "dip" angle is finite, 
therefore the QO^f signal originates from a smal] segment of tlie PPT plume and 
a segment of the ^^0LTRAP wall. The insert in Figure 1 indicates the kind of 
data expected from such a measurement. The greatest backflow would be expected 
at a dip angle of 0 ° - decreasing to lower values in the downstream direction. 

Once the view angle begins to intercept the plume boundary-wall intersection, 
the wall backscatter will begin to dominate the signal. By considering only 
the data for small dip angles, the net total plume backflow can be calculated 
by sunning the data over the dip angle distribution. Using a collimator with 
a conical aperture as shown in Figure 2 would require observations over various 



Figure 1. Backflow Measuring Technique 





dip angles both parallel and perpendicular to the thruster axis in order to 
observe the entire pluiw? volume. To avoid the experimental con^ilexity 
associated with such an aperture, a slotted aperture was chosen for the collimator 
design, as also shown in Figure 2. 

Building upon the previous Phase I efforts, the Phase II investigations 
of the PPT primary plume and plune backflow have been completed, and are 
detailed in this report. A discussion of the primary plume measurements and 
their impact on the plume backflow testing is followed by a description of the 
design of the Q(31 backflow ttst apparatus. After these sections, tlie 
collimated QQ4 measurements taken in the MOLSINK facility are described and used 
in a detailed analysis of the PPT plune backflow level and distribution. Finally, 
an additional section is includSv’ which describes a small follow-on measurement 
of the plume backflow from the PPT with a radically modified nozzle geometry, to 
determine the backflow sensitivity to nozzle design. 
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2.0 PRIMARY PLIM; STliDIRS 


The source of the backflow from the pulsed plasma tJirustor is the primary 
plume downstream of tlie thruster nozzle. The backflow mass flux magnitude 
and its distribucion away from and around the thruster axis both d,-pend on the 
[iriiruiry plume mass flux, velocity, and chemical structure. In order to develop 
an accurate picture of tlie PPT backflow, and, more practically, to assist in the 
design of an acceptable backflow measurement teclmique, a basic understanding of 
these characteristics of the PPT primar>' pliime is necessar>^ To develop this 
understanding, various tests have been carried out and are described in tlie 
following paragraphs. Tlie major fraction of these tests was performed in the 
ShP vacuum facility duo to its greater operating convenience and lower cost . 

The remaining tests require a thruster environment more similar to actual 
space flight conditions, and, hence, were performed in the KDLSINK facility. 

2.1 Mylar Target Deposition The exit orifice of the PPT exhaust nozzle is a 
11. f) by 16.5 an rectangle with a re.sulting aspect ratio of 1.4. This azimutlially 
nons>'minetric shape iniplies that the exhaust plume, and hence the plume backflow, 
may also lie non-axi symmetric around tlie thruster axis. In addition, any plume- 
wall backscatter would bo non-axisymmetric. This would require any primary or 
backflow plume measurement to bo made at various azimuthal as well as radial and 
axial locations, and, thus, would considerably complicate the experimental testing. 
Recent observations at the Fairchild Republic Co.^^^ indicate that 40 cm 
downstream of the nozzle exit plane, the PPT plume is elliptical in cross- 
section with its major axis parallel to the nozzle longer dimension, and with 
an aspect ratio of only 1.2. These results suggest that at greater domstream 
axial locations, the plume may approach azimuthal s>'mmetry, 

IVlien installed in the MOLSINK facility for the plume backflow measurements, 
the PPT was placed approximately in the center of the enclosed volume so to 
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provide maximum themal isolation frotn the MOLSINK walls. In this position, 
the thruster plume will collide with the MOLSINK walls approximately 80 cm 
downstream of the PFf nozzle exit plane (assuming a plune expansion angle of 
about 40”^^^). To check the primary plume symmetry, and hence the plume-wall 
backscatter s>Tnnetry at this axial location, a technique developed for plune 

studies cf the 8 cm ion bombardment thruster was used. A 1.2 meter square 

0 

sheet of 1 mil thick Mylar coated with a 700 A thick layer of aluminum was 
placed 76 cm downstream of the PPT nozzle exit plane, on a frame siqjporting its 
top and bottom edges. A 38 cm diameter hole was cut in the center of this sheet 
to permit the central core of the PPT plume to escape without damaging the 
fragile f»fylar. A photograph of this aluminized iV^lar target is shown in figure 3, 
as installed in the SEP vacuum facility prior to the test. 

During the test, the target could be seen to oscillate after each thruster 
discharge due to the plume impingement. The target was exposed to approximately 
12,000 discharge pulses over a three day period, and then removed from the tank 
for analysis. A photograph of the target after the test is shown in Figure 4. 

A series of concentric rings can be observed in this black and white reproduction 
which, in fact, are multicolored, in a manner similar to the bands of light seen 
in quarter-wave diffraction plates. Such plates consist of a, highly reflective 
surface, such as aluminum, covered with a layer of transparent material with 
a thickness equal to an odd multiple of a quarter of a wavelength of the 
absorbed light. Thus, the presence of these concentric rings on the aluminized 
target indicates that material has been deposited. 

Tlie primary reason for performing the plume target test was to determine 
the azimuthal symmetry of the PPT plume, 70-80 cm downstream of the nozzle. The 
shape of the concentric rings seen on the plume target provides an accurate 
measure of this symmetry. The center of these rings is displaced upward 
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approximately 9 an from the center of the plume target, indicating a slight 
misalignmoit in either the plume discharge or the thruster mounting. The 
concentric rings are circular, except for flat spots at the 10 o’clock and 
2 o’clock positions. The top of the target (12 o’clock) corresponds to the 
location of the cathode elec.trode and the spark plug trigger; hence these 
flat spots may be associated with the difference between the cathode and the 
anode discharge physics. In any case, these flat spots represent a deviation 
from the circular ring mean radius of less than 5%; thus for practical 
purposes the primar)’’ plume can be taken to be azimutlially uniform. A further 
test to check the backflow plume axisymmetry is described in section 4.2 and 
confirms that the backflow is also azimuthally uniform. 

In order to obtain more detailed information on the plume deposit, three 
different methods of measuring the actual thickness of the deposit on the 
plume target were attempted. The first method, utilizing a laser ellipsometer, 
failed due to the flexibility of the plume target M[ylar substrate. Tliis flex 
ibility prevented the target from resting evenly on the sensing platform, which 
led to large inaccuracies. The second method utilized a Sloan DEKTAK Surface 
Profile Measuring System which senses the position of a scribe as it moves along 
the sample surface. This method was also unable to measure the absolute 
thickness of the deposition, again because of the flexibility of the Mylar sub- 
strate; however it did indicate tliat the existing micropores in the aluminum 
layer were smoothed over towards the target center, llris could happen either 
by the deposit filling the micropores or by the plume eroding the surrounding 
aluminum. Finally, the transmittivity of the target was measured to see if any 
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qualitative evidence could be found to determine the deposit thickness. The 
transmission coefficient was found to drop from 2.7% at the outer edge of the 
target to 2.0% at the edge of the hole in the center of the target. This 
change in transmitt ivity is relatively small and indicates that either the 
deposit thickness increases towards the target center or that the target 
surface features cliange with decreasing radius to increase the reflectivity. 

In the light of the smoothed over micropores found with the surface profile 
measurement, it is more probable that the target surface features changed. 

None of the previous methods were able to determine the difference between 
erosion or deposition on the target, and hence, provide little improvement 
over the previous quarter wave plate analysis, in the understanding of the 
plume target results. 

rs 10 ii") 

2.2 Double QCM Probe Mass Flux Previous investigators '•*»-' have studied 
PPT plumes using such various diagnostics as Langmuir and B-field probes, 
calorimetric disks, high speed photography, glass capture caps, Faraday cups, 
microwave interferometry and single QCMs. When combined with the known per- 
formance of the PPT, these studies indicate that a significant fraction of the 
plume mass flux consists of low energy (probably neutral) particles. Tlius, 
in order to measure the radial distribution of the PFf primary plume mass flux, 
a method is needed \diich is sensitive to both the neutral and ionized components 
of the plume. In addition, the method must provide adequate spatial resolution 
and an in situ, real time output to minhnize error. 

To .'latisfy the above requirements, a Q01 measurement would seem to provide 
an adequate solution; however, as previously mentioned, earlier attempts 
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to determine the radial distribution of the plume mass flux with a single QCM 
have met with little success because the plume erodes the QCM collecting surface 
rather than depositing on it. In order to alleviate this difficulty and still 
maintain the advantages of a QCM measurement, a special double QCM probe was 
developed, and is sketched in Figure 5. It consists of a shielded container 
with an aperture designed to direct the incoming mass flux to QCM 1, placed 
at an angle of 45° with respect to the incoming axis. The material which 
reflects or erodes from CJCM 1 is partially captured by QCM 2, placed normal to 
the incoming axis and on the optical path from QCM 1. 

The net signal output, S^, of QCM 1, facing the PPT plume is equal to 
tile rate of material deposited on the costal sensor. This rate is equal to 
the axial component of the local plume mass flux, ih, less the amounts 
reflected and ablated, from the QCM surface: 

Sj = ill cos 45° - m^ (1) 

The cos 45° is required to correct for the angle between the incoming axis and 
the QCM collecting surface normal. The reflected and ablated mass flux leaves 
the surface of QCM 1 in some unknown distribution about the QCM surface normal. 
Some fraction, K, of this mass flux impinges on and is collected by QC3^1 2. 

Thus, the signal output, S2, of QCM 2 is: 

S 2 = (2) 

The fraction, K, not only accounts for the fraction of reflected and ablated 
material from QCM 1 that impinges on QQ4 2, but also for that fraction of 
impinging material that actually sticks to the collecting surface of QCM 2, 
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Figure 5. Double QCM Probe Schcanatic 







rather than ablating and reflecting from it. 

The axial plume mass flux, A, can be found by eliminating the reflected 
and ablated mass flux, from Equations (1) and (2): 


A 



( 3 ) 


The values of and S 2 can be measured locally to give a local value of 
A, provided the constant, K, is known. This constant can be found through a 
caloulation of the total mass flow rate, rfUp, from the thruster discharge as 
follows . 

The total mass flux over the entire plume cross section is given by: 



(43 


where dA is a differential cross-section area element, normal to the thruster 
axis. Substituting into equation (3): 



Under the assumption that tlie constant, K, is independent of radius, 
equation (5) can be solved for K: 


K = Aj. cos 45° ~ 

7 


dA 


( 6 ) 


rS2dA 


'Hie total mass flow rate, rfij,, is knoMi to be 1.56 mg/pulse; thus 
measurements of and $2 versus radius can be used to experimentally 
evaluate the constant, K. Once known, K can be substituted along with 



local values of Sj^ and S 2 into equation (3) to give the local axial plume mass 
flux, ih. 

The double QCNl probe used in the testing is shown in Figure 6 with and 
without its cover plate. The quartz blank of QCM 1 with both its collecting 
and reference electrodes can be seen. QCW 2 has a separate cover plate to 
insure that the mass collected on QCW 2 is only from that eroded and reflected 
from the QCM 1 collecting electrode, and not from any spurious internal 
scattering. The circular shaft shown in the photographs is a mounting fixture 
used for assembly that simulates the required liquid nitrogen (LN2^ cooling 
line. IVhen installed in the SEP vacuum facility, the IJ ^2 cooling line 
consisted of a 1.5 inch diameter flexible stainless steel tube fed through 
the top of the vacuum tank and down to the double QCM probe where its end was 
plugged. This cooling line was filled with LN 2 outside the vacuum tank, 
while gravity acted to keep the LN 2 probe, 

The 1.5 inch diameter was necessary to prevent a vapor lock from forming 
and preventing the LN 2 from reaching the double QCM probe.. 

The double QCM probe was mounted on a movable support that was capable 
of sweeping the probe radially outward fron the thruster axis to a radius of 
roughly 75 cm. A complicating factor in the construction of this support was 
the requirement that the probe remain tied to the LN 2 feed line. Although 
made of flexible steel tubing, at LN 2 temperatures this line is relatively 
stiff, requiring that tlie support he sturdy enough to move the probe agaiiist 
the drag of this LN 2 iii^e. Two design options were considered: 1) a support 
rack which pivots about a point in the PPT nozzle exit plane at the thruster 



raSS" : ' ■''■•■Ate 



Hi ’ 








axis, such that the double QCM probe vas always at a constant distance from 
this point; and 2) a support rack which moves radially along a straight line 
perjjendicular to the thruster axis, such that the probe remained a constant 
distance downstream of the PPT nozzle exit plane. Tlie first option would have 
reqiiiral the construction of a strong, curved track and a complex motor assembly 
to move the support along this track. In addition, maintaining the alignment 
of the probe would have been difficult. Because the second option is meclianically 
more simple, it was chosen for the support design. 

Figure 7 shows a photograph of tlie final installation in the SEP vacuum 
facility. The probe is mounted in a support which slides on Teflon bearings 
along two parallel stainless steel tubes. A cable -chain drive system is 
connected to the probe support around two pulleys, seen at each end of the two 
parallel tubes. This system is driven by an electric motor at the base of the 
diamond shaped structural frame. The LN 2 feed line can be seen curving from 
the probe up to the top of the vacuum tank. With this system , the douljle QCM 
probe can be positioned anyvdiere between 30 cm to the left of the thruster axis to 
75 cm to tlie right. The probe cover plate was positioned 74 cm downstream of 
the PPT nozzle exit plane, approximately at the same location as the aluminized 
Mylar target discussed in the previous section. 

Using the double QO^ probe, measurements were made over a period of several 
days with the thruster firing once every 17 seconds. The SEP facility walls 
were maintained at LN 2 temperature to minimize the wall -plume backscatter, Tlie 
data taken during this test were reduced to the mass flux values, and S 2 , 
for both (y’M.s in the probe. These are shown in Figure 8, plotted versus radius 





Figure 8. Double QCM Probe Data 
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measured from the PPT plume axis. As expected, the signal from Q04 1 indicates 
its collecting surface is eroding for probe positions out to a radius of about 
45 cm. The signal from QCM 2 is consistent with this erosion in the sense that 
it is largest wlien the eroded mass flux is largest, and hence when the great- 
est amount of material is available for collection. 

'fhe two data points shown for each QCM at a radius of 26 an were taken 
on opposite sides of the thruster axis, and thus provide a measure of the 
symmetry of the plume and tlie accuracy of the probe analysis. The two data 
points of the QCM 1 signal are virtually identical, indicating that the incoming 
^lurne mass flux is the same on opposite sides of the PPT axis, and, hence, is 
ax i symmetric. The two data points of QCM 2 differ by about a factor of 1,6, 
indicating that the fraction, K, of scattered and eroded material collected 
by QCM 2 from QCM 1 can vary by as much as 60%. 

Using the data of Figure 8, Equation (6) was used to calculate tlio constant, 
K. lliis value was found to be 0,014, and indicates that only 1.4% of the 
material reflected and eroded from QCM 1 is collected by QCl^l 2. The error 
induced in the calculated mass flux profile due to the variation in K can be 
seen in Figure 9 which shows this calculated profile versus radius. The two 
data points at a radius of 26 cm indicate that the error in this mass flux 
due to the variation in K can also be as much as 60%. With this implicit 
error in mind, the profile of Figure 9 still provides a reasonable measure 
of the mass flux distribution over the plume radial extent. This profile 
approximates a Gaussian shape, witii a half width at half maximum of 28 cm, 
corresponding to an enclosed half angle from the plume axis of 21°. 
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At a radius of 60 cm (wliich corresponds to an enclosed lialf angle of 40°) the 
plume flux is less than 1C% of its centerline value. 

2,3 High Speed Plume Photography Earlier studies of lower energy pulsed 
plasma thruster discharges have utilized high speed photography to study 
the evolution of the primary plume over the total discharge time. These 
studies have provided valuable information about the PPT primary plume 
velocity and formation, which indicates that the plume is not a simple 
expansion of a homogeneous plasma. In order to investigate the properties of 
the 1 mlb PPT plune under study in this report, and to substantiate some of 
the previous experimental results, a high speed photography study of the PPT 
primary plume was carried out. In particular, a measure of the plume velocity 
was desired in order to determine the time at which the plume collides with 
the vacuum tank wall and begins to backscatter. 

This study was carried out with the thruster installed in the SEP vacuum 
facility in the same location as for the previous >fy^lar target and double 
QCM probe tests. An observation port in the side of the vacuum chamber was 
used to view the plume. This port has a glass window with its center displaced 
46 cm downstream of the PPT nozzle exit plane. The experimental set-up is 
sketched in Figure 10, and shows the trigger delay generator (TRW model 46A) 
which served to trigger the camera shutter after a set delay from the begin- 
ning of the PPT disdiarge pulse. The camera is a TRW image converter camera 
model ID and used a microsecond framing plug-in unit Model 4B to control 
the number and length of exposures per discharge. Tliis particular framing 
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unit allowed three separate exposures to be made during each PFT disdia>'ge. 

Tlie duration of each exposure was of order 0.2 psec and the separation 
between each exj^osure was controllable between 0.5 and 20.0 ysec. To insure 
that the axis of tlie camera was perpendicular to the PPT plume axis, and to 
determine the exact image demagnification, a ruler was suspended from the top 
of the vacuum tank, along the thruster axis. The camera was aimed at a point 
on this ruler corresponding to a right angle between the thruster and camera 
axes and a photograph of the scale was taken. Tlie demagnification was tlien 
calculated and found to be 13.9:1. 

Figure 11 shows a sample series of photographs taken of the PPT plinne. 
flacli individual frame has three grid lines superimposed on the actual plume 
exjiosure. The overall series shows the time history of the plume as it 
passes the observation area, with time increasing from the bottom of the 
figure. Each group of three pictures was taken during one discharge pulse. 

In order to compare the plume behavior from pulse to pulse, each separate 
series of photographs overlaps the time of the accompanying series. The 
photographs in Figure 11 span the time from when the plume first appears in 
the viewing region, about 19 psec after discharge initiation, to when the 
plume luminosity decays to where it is no longer visible, about 37 psec 
after the discharge initiation. As observed, the plume is not homogeneous, 
but in fact contains locallized regions of higli luminosity, and thus 
presumably high density. Over fifty photographs similar to those in Figure 11 
were taken aiid all show some degree of nonunifonnity in the luminosity pattern. 
The nonuniforinities in the plune plasma must average out over many discliarge 
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pulses, in order to give the relatively uniform experimental results found 
with the Mylar target and double QCM probe tests. The plume backflow is 
expected to behave in a similar manner, since it originates in the primary 
plume. 

The velocity of the luminous disdiarge plasma seen in the photographs 
of Figure 11 can be detennined by calculating the time -of -flight distance 
between one pliotograph and another, separated by a known time interval. 

Using this method, the velocity of the plasma front as seen in the earliest 
group of pictures in Figure 11 is 30 + 5.5 kra/sec. The velocity of the 
bright luminous region, seen in the middle group of Figure 11, is 23 ± 5.5 km/sec. 
The error bars in tliese velocities represent standard deviations calculated 
by combining the measured velocities of several series of photographs similar 
to those in Figure 1 1 . The discrepancy between the velocities of the plume 
leading edge and the lujiiinous region inside the plume may be due to several 
reasons, including different magnetic force acceleration patterns, and 
different local acoustic flow properties. 

The plume average exhaust velocity is known to be 17 km,/sec, which is 
significantly smaller than the tut) measured velocities. Tliis suggests 
that the luminous portion of the plume consists of only a fraction of 
the total discharge mass and that the remainder is moving at a velocity 
lower than the pltmie average. This t)’|)e of behavior has been seen before 
in micro]X)und thrusters , where the luninous portion of the plume was 
found to be tlie ionized fraction of the plume. Measurements of the ion 
velocity of the millipound thruster plume were made at Fairduld using 
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a Langmuir probe and sliow good agreement with the present photographic 
measurements. 

In the M3LSINK facility, where the plume backflow measurements were made, 
the thruster was installed with its exit plane perpendicular to the major axis 
of the ellipsoidal tank, roughly l.S m from the tank end. Thus, the plume 
leading edge will reach the tank wall roughly 50 ysec after the PPT discharge 
initiation. Since this portion of the plume is the most energetic, it would 
be expected that most of the wall backscatter would originate with this 
portion. Assuming, conservatively, that the backscatter velocity is equal to 
the incoming plume velocity, the backscatter will reach the area of the thruster 
nozzle roughly 100 psec after the discharge initiation. 

Hie energetic portion of the PPT plume is ionized and thus is confined to 
within the magnetic field pattern of the PPT discharge. The lower energy, 
slower moving portion is not confined by this field and thus, is free to flow 
radially outward and axially upstream more easily. Hie PPT plume backflow is 
expected to be primarily composed of material from this lower energy portion. 
Although the velocity of this portion of the plume is unknown, it must be 
less than the plume average velocity and is probably close to its sonic 
velocity. Assuming the temperature of tliis material is less than 10,000 K, 
an upper estimate of the velocity can be found from the definition of the 
sonic velocity, Cg! 

Cg = V Y RT 
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where y is the ratio of specific heats and is taken to be 5/3, and R is the 
ideal gas constant divided by the average plume molecular weight of 16.7 amu^®\ 
ITie result gives a value of roughly 3000 m/sec, or about one-tenth of the high 
speed plime portion. Using the estimated wall backscatter return time of 
100 ysec, the low speed backflow will travel only about 30 cm from the thruster 
axis before the wall backscatter overtakes it, thus any backflow measurement 
in the MOLSINK facility must be designed to correct for an almost simultaneous 
wall backscatter. 

In addition to the preceding observations of the PPT plume on the thruster 
axis, the high-speed camera was also used to determine the expansion angle of 
the luminous portion of the plume. Figure 12 shows a schematic of the camera 
geometry used to make this measurement. At the camera axial position of 46 cm, 
it was tilted up from its original position, aligned perpendicular to the 
thruster axis, to a position where the edge of the plune was centered on the 
photograph. As shown in part (a) of Figure 12, the angular displacement was 
found to be 17”. Using ttos angle and the known distance from the 
camera to the thruster axis (132 cm) , the plume radius at this axial location 
can be determined. As shown in part (b) of Figure 12, the plume expansion 
angle, e, can then be calculated and I'oiuid to be rough!. VO”, which is 
in agreement witli t!ie previous measurements. From the photographs taken of 
the PP! plume edge, the velocity of tlie plume at tliis location was found to 
be 26 + 5 km/sec. Tliis velocity is approximately equal to the measured center- 
line plume velocity; thus the plume velocity radial profile is essentially flat 
out to the plume edge. This type of radially unifoim velocity profile is 
similar to those found in other t)'pes of plasma thrusters. 



I 


j 



(a) END VIEW 



(b) SIDE VIEW 


Figure 12. PPT Plume Camera Geometry 
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2.4 Plume Composition Analysis The results of the previous studies have 
indicated that the plume material in the central part of the primary plune 
is energetic enough to erode the vacuun tank wall surfaces. In addition, it 
is loiown that the major amount of plume -wall backscatter is from this central 
part of the plume. These two facts combine to suggest that the backscattered 
material from the tank walls may be composed of ablated material from the wall, 
and hence, recognizably different from the backflow from the plune itself. In 
order to determine if this possibility is indeed true, a series of tests was 
run to determine the plune composition using visible light spectroscopy and carbon 
disk analysis. 

A 0.5 m Jarrel-Ash grating spectrometer was set up in the same location 
as the high speed camera shown in Figure 10, to observe the optical radiation 
from the PPT exhaust plume perpendicular to the plume axis. The wavelength 
range from about 2000 K to 6000 K was covered using Royal Pan Film, with a 
n»ercury vapor lamp for a comparison spectrum. Exposures were varied frcxn 
15 to 45 discharge pulses to provide adequate resolution. The results indicate 
tliat spectrfi lines can be found over the entire range. Analysis indicates 
that much of this radiation is due to singly ionized fluorine with some 
contributions from ionized carbon. No lines from neutral carbon or fluorine, 
or from any other specie, were found. Tliis indicates that the energy of the 
plume neutral specie component is such tliat the radiation from tliis component 
is negligible with respect to that from the ionized specie. It also suggests 
that little reccanbination of the charged particle plasma is occurring upstream 
of the observed plume region. Since no spectral lines from any specie from 
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the vacuvan tank wall were observed, it appears that the wall backscatter 
cannot be resolved from the plume backflow using this spectroscopic method. 

An alternate method of resolving the possible differences between the 
plume backflow and wall backscatter involves the analysis of the material 
deposited on sanpling surfaces exposed to the PPT plume. A scanning electron 
microscope (SHM) was used to observe the morphology of the surface deposits 
on these san^jling surfaces and energy-dispersive X-ray analysis was used 
to identify the atomic species. This type of X-ray analysis is unable to 
differentiate between species of atomic numbers lowr than 9; hence only 
the fluorine in tlie Teflon propellant can be identified, idiile the carbon 
will remain transparent. In addition to possibly differentiating between back- 
scatter and backflow, scanning electron microscopy is also useful in diecking 
the uniformity of the deposits on the collecting surfaces. This uniformity 
strongly a ffects the calibration constant for the Q04s used in these plume 
studies. For the previous work, including tliat of the Phase I effort, this 
deposit was assumed to be uniform and the calibration constant was calculated 
accordingly. This assunqition can be checked using the aforesaid method. These 
analyses, including both the SEM photography and the X-ray energy dispersive 
spectroscopy, were performed by Dr. Raymond L. Chuan of tlie Brunswick Corporation. 

Initially the analysis was carried out on the deposit on the collecting 
surface of a QCNI used in the off-axis skimmer of the Phase I testing. 

It was hoped that an examination of tliis QC^l would provide sane evidence 
of the species backscattered from the tank wall, since only this backscatter 
could have reached tlie QC^l. An SEM pliotograph of a portion of this QCM 
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doilltecting surface shM In Fi|;ure 13 (a) . As evMiglll%| th% deposit on the 
■^(3^ platinim electrode siirfS^e to^nsists of relatiVfely Wm isolated particles . 
Particle A is shown at a gTeatdr magnification in 13 (b) , where it 

aplpcars to be amorphous with an overall sizfe of roiighl^ 80 microns. X-ray 
spectra of this particle indicate that it is ctW 5 >osbd primarily of alumihiim 
and hence, is probably a sputtered particle from the aluminum in the vacuum 
tank walls or QQl body. Thfe remaining particles were each examined and 
found to have no resolvable X-ray spectra, indicating lhat they are composed 
of low atomic number elements (less than 9) , Which may or may not be from 
the PPT plune. 

Analysis of the previous QCM surface was ambiguous, sinca its collecting 
surface was shielded from the PPT plune by the skimmer wall and since little 
measurable mass was deposited on its surface. In addition, the X-ray spectra 
of the few particles on this QC^l surface were masked by the spectra of the 
silicon from the quartz crystal and the platinum from the actual collecting 
electrode. To remedy these problems, three carbon disks approximately 
1.0 cm in diameter were installed in the MDLSINK facility to be exposed to 
the PPT discharge. Carbon disks were used because carbon has an atomic number 
of less than 9, and hence is transparent to the X-ray spectroscopy used. 

Oadi disk was placed in the bottom of a 2.0 cm by 3.0 cm box, roughly 1.0 cm deep, 
which acted as a relatively open collimator to control the region viewed by the 
carbon disk surface. These boxes were either pure aluminum foil or commercial 
pot metal depending on the particular carbon disk. Two of the carbon disks 
were placed side-by-side on the dbwhstream edge of the thruster aluminium 
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enclosure, approximately IS cm off the thruster axis and 5 cm upstream 
6f the exit plane. They iwre set to face directly downstream toward the 
M3LSINK wall area where the backscatter is greatest. One of these disks 
was mounted in a pure aluminum foil box and the other in a pot metal box. 

The third disk was attached to a bracket on the MDLSINK wall, about 45° 
off the thruster axis. This disk was set to look directly into the thruster 
discharge chamber and was in a pot metal box. These disks were exposed to 
approximately 110,000 discharge pulses and then removed for analysis. 

The analysis of the two carbon disks mounted together on the thruster 
enclbsure indicates that there is a quantitative difference between the 
deposit collected on the disk in the pure aluminum box and the deposit on 
the disk in the commercial pot metal box. This difference indicates that 
some material was actually eroded from the boxes containing the carbon disks. 
Thus, any elements vdiich are contained in the box material and are seen in 
the analysis may not necessarily originate in the PPT plume or tank wall back- 
scatter. These elemaits include aluminum from all the boxes and iron, lead, 
zinc, and other trace elements from the pot metal boxes. 

The carbon disk from the pure aluminum box mounted on the thruster enclosure 
and facing downstream was analyzed to determine the deposit morphology and 
elemental composition. Two sample SEM photographs of the carbon disk surface 
features are shown in Figure 14. Figure 14 (a) is a low magnification view 
idiich shows a number of particles adhering to the surface. Figure 14 (b) 
is an enlarged view of particle A showing it to be an amorphous lump of 
material seemingly composed of many 
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(b) PARTICLE A 


Figure 14. SBl Photographs of Carbon Disk Deposits 
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small spherical particles. X-ray spectroscopic analysis of this particle 
shows a strong fluorine line which is probably in chemical combination 
with carbon, making this particle similar to Teflon. Figure 15 shows two more 
enlarged views of particles observed on this carbon disk. X-ray analysis 
indicates that these particles also show Teflon; however, their morphology 
is fundiunentally different from that of Figure 14 (b). Both of these particles, 
and in fact almost all the observed particles, appear to be cracked from 
a uniform layer built up on the carbon disk substrate. Tliis cracking may 
possibly be due to thermal stresses induced during the warming of these disks 
to room temperature after having been maintained at low tenperatures during 
the experiment. 

Although the carbon disk mounted on the MOLSINK wall was installed in a 
pot metal box, and therefore subject to considerable contamination, the 
analysis of the deposit on this disk indicated several interesting points. 

To the naked eye, this disk was well covered with a material showing a 
velvety purple color similar to the deposits seen on the MDLSINK walls 
around the lower door. The surface density of particles on this disk was 
higher than that of the other disks; however, the analyses of these particles 
indicate that they are composed of materials from the pot metal box, and tlius 
may not be due to the PPT discharge plume. As in the previous disk, this disk 
also showed the ubiquitous presence of fluorine over the entire surface, 
suggesting that the teflon propellant deposited in a unifoim layer. In addition, 
copper, presumably from the PPT electrodes, was also seen to be uniformly 
distributed over the disk surface. This indicates that it is possible to 
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Figure 15. SEN! Photographs of Surface Particle Cracking 
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qualitatively recognize the PPT discharge erosion products in the discharge 
and, hence, possibly determine the actual erosion rate. 

Hiis analysis of the carbon disks exposed to the PPT discharge indicates 
that there is no easily recognizable difference between the material rebound- 
ing from the M3LTRAP wall and the material flowing directly from the thruster 
discharge chamber. Thus, it would be impossible to distinguish between the 
PPT plane backflow and the plume -wall backscatter using this method. Further 
analysis indicates that the major part of the deposition on these carbon 
disks is in a uniform layer with only a few particles at isolated points. 

Tfiis type of deposition indicates that the earlier assunption of uniform 
deposition on QQ4 collecting electrode surfaces is accurate; hence the use 
of this fact in calculating the QCM calibration constant is justifiable. 

2 . 5 Sumnary The previous studies of the PPT plane have provided valuable 
insight which can be applied to the design of an appropriate plume backflow 
measuring system. The measured axisymmetry of the plane downs treajn of 75 cm 
from the nozzle confirms that the plume-wall backscatter is axisymmetric, as 
was indicated in the Phase I studies. Furthermore, this axisymmetry suggests 
that the PPT plume backflow may also be axisymmetric; hence the experimental 
program to measure this backflow need not include an extensive study of the 
azimuthal variation in this backflow. The radial mass flux measured using 
the double QCM probe indicates that virtually all of the primary plume is 
confined to a 40° half-angle conical expaniion. This measure of the plume 
boundaries and the estimate of the flux density within this plume will be 
useful in determining the regions to be observed in order to measure the plume 
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backflow. Finally, the plume velocity and composition analysis indicate that 
the plime backflow cannot be easily differentiated from the plume-wall 
backscatter, by either appropriate sensor timing or elemental analysis. Thus, 
to measure the plume backflow, a method must be developed which differentiates 
between the backflow and the backscatter in sane other fashion, such as the 
method described in Section 1, using collimated QCMs, 
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3.0 BACKFLOW MEASUWM3m’ DEVELOPMENT AND TEST HISTORY 


Based on the results of the previous section, it appears that the 
conceptual method of measuring the PPT backflow using collimated QCMs is 
feasible. This method requires the design and assembly of an experimental 
ypparatus capable of supporting an array of collimated QCMs, moving these QOIs 
to various radial and axial locations in the plume, and varying the dip angle 
between the PPT nozzle exit plane and tlie collimator axis. In addition, eaclv 
QCM must be cooled to LN 2 temperatures and have its temperature regulation, 
power, and output signal leads connected to tlie appropriate systems outside of 
the MOLSINK tank. The first half of this section will discuss the overall 
design of this experimental apparatus. Included in this discussion will be 
a description of the preliminary testing of various collimator designs carried 
out in the SEP vacuum facility. 

After tlie assembly of the test apparatus, it was installed in the 
MILS INK facility and used to gather the necessary experimental data required 
for tlie determination of the PPT plume backflow. The test history and the 
reduction of tlie test data from QCM beat frequency shifts to mass flux values 
will be discussed in the latter half of this section. 

3.1 Collimator Design and Testing The conceptual technique for measuring 
the PPT plume backflow was discussed in Section 1.0. Referring to Figure 1 of 
that section , the collimated QCM signal consists of contributions from the PPT 
plume and the MOLTRAP wall area within the collimator observation region. For 
small dip angles, the plume-wall backscatter can be ignored; thus the collimated 
Q01 signal ivould he duo to the hackflow from the observed voliune of tJic PFl' pliune. 
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By using a collimated QCM which observes a plume slice of enclosed angle, 

Ant , ^d observing these slices from a dip angle of zero out to the maximum 
yalue for negligible wall backscatter in steps of Aa , the partial plume 
backflow from this total volime might be found by simply siinming the 
measured signals. Unfortunately, physital limits of the collimator design 
prevent this simple procedure from giving accurate results, as can be seen 
by the following. 

I/)oking in a direction perpendicular to the QO^I collimator axis, the 
geometry is as shown in Figure 16. The region observed by the QCM collecting 
electrode can be divided into two subregions: the illuminate and the penumbra. 
Any point source of backflow in the illuninato will see the entire electrode 
surface; hence the measured QQM signal will be directly proportional to the 
electrode area. Any point source in the penumbra will see only a fraction of 
the electrode surface (due to shadowing by the collimator aperture \ and hence, 
will depend on the electrode area in a more complex fashion. This partial 
shadowing must be corrected for in the summation of the QCM signals at different 
dip angles in order to insure that all the backflow over a given range of dip 
angles is measured. This correction is analytically very complex, as will be 
seen in Section 4, and hence, it is desirable to design the QCM collimator such 
that the penumbra and the associated correction to the measured QCM signal is 
small. As will be seen by the following analysis, this cannot be done under 
the existing experimental constraints. 

Referring to Figure 17, the edge of the penumbra is at an angle, 0 , with 
respect to the collimator axis, and the edge of the illuminate is at an angle 


- 46 - 







Figure 17. QQi Collimator Aperture Geometry 
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with respect to this axis. With the QCM electrode width, q, and the 
electrode-aperture distance, s, the penumbra angle, 3, is: 


tan 


-1 


s 


+ tan 


Aa 

T- 




From Figure 1, it can be seen that the aperture angle , Aa , should be less 
than 10° to insure adequate spatial resolution. Table I shows the values 
of calculated for various Aa and q/s. The QCll electrode width, q, is 
approximately 0.8 cm; hence the tabulated values of q/s cover a range of elec- 
trode-aperture distances, s, from 8 to SO cm. 


TABLE I. QOl Collimator Penumbra Angle (degrees) 


Aa 

(degrees) 

q/s 

0.1 

0.05 

0.01 

0 

S.71 

2.86 

0.573 

5 

8.17 

5. 35 

3.07 

10 

10.6 

7.82 

5.59 

20 

IS. 4 

12.8 

10.6 


f^ractical considei'ations of the available space in tlie ^DLSINK tanlc 

dictate that s should be no larger than about 10 cm; hence the 

ratio q/s is restricted to values greater than rotjghly 0.1. According to 

Table I, this indicates that the collimator penumbra angle will be equal to or 

larger than the aperture angle; hence the correction tc the collimated QCM 

signal due to the penumbra must be large. 
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The choice of the collimator aperture angle depends on a con^rcMise 
between the good spatial resolution of small apertures and the magnitude of 
the QCM signal which decreases with small apertures. The QQf signal not only 
depends on aperture size but also on exposure time, dip angle and location. 

In general, the signal decreases with increasing dip angle (out to the angle 
where the wall backscatter begins to increase) and with increasing distance 
from the thruster axis. To determine the magnitude of the QCM signal* and 
hence, aid in the final dioice of collimator aperture angle, a series of tests 
was run in the SEP facility with two types of QCM collimators at various 
locations. The collimators are identical except for the aperture angle, vdiich 
is 20° for some collimators, and 0° for the remainder. A perspective 
view of the 0° aperture angle collimator is shown in Figure 18. The 
QCM crystal is ex^iosed to the main collimator through a 0.8 cm square hole 
cut in the center of the collimator backplate. The front face of the 
collimator is circular so tliat the apertui'e angle remains constant over the entire 
width of the collimator slice. This width is set at 100° to include 
the entire width of the PFT plume and yet not over expose the regions outside 
of the plume idiich would contribute to the observed wall-plume backscatter. 

The radius of the curve front is 8.0 cm, exactly 10 times the QCM electrode 
width; hence the penumbra angle for the 20° collimator is 16° 
and for the 0° collimator is 5.7° (see Table I). 

Six collimated ()CMs were mounted in the thruster nozzle exit plane in 
a rectangular array as shown in Figure 19. The three rows of two QCMs each 
were placed 48 cm, 63 cm, and 78 cm from the thruster axis, respectively. 


COLLIMATOR SLIT 



Figure 18, QCM Collimator Design 
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Figure 19 . QCM Collimator SEP Facility Test Arrav 
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All the collimators were set to view along a line cutting the thruster axis 
at a dip angle of 60"^. The two QCMs at each radial location were expected 
to be insensitive to their slightly different azimuthal positions due to the 
measured plume symmetry about the thruster axis. As will be seen, this assuiiqj- 
tion is acceptable within the error of the measurements. In Figure 19, QGM 
numbers 3 and 5 had the narrow 0° collimators wliile QCM numbers 1,2,4 
and 6 had the collimators. 

In the SEP facility, the PPT thruster was fired downstream towards a large 
LN 2 cooled steel target at the end of the tank. This target can be rotated 
about a horizontal axis, perpendicular to the thruster axis, and thus, was used 
to vary the tank backscatter characteristics during the test. In addition, the 
Q04 array was in place during the previously discussed plume Mylar target test, 
and data was taken during tiiis test. The results of the QCM collimator 
testing arc shoivn in Table II for the six QQ^ and the various test conditions. 

The positions and dip angle of the collimated Qd'ls were chosen to minimize 
the plume-wall backscatter effects on the observed data. In the SEP facility 
it is clearly impossible to eliminate the backscatter, and the data of Table II 
can be used to determine the magnitude of the backscatter effect. The data 
taken with the pliuiic Mylar target in place is generally about a factor of two 
larger than the data taken without the target, indicating that a large fraction 
of the plume is being backscattered by this target. This backscatter increase 
is essentially independent of QCM collimator aperture angle, but does seem to 
increase with increasing QCM radial position. The downstream location of the 
Mylar target ras chosen to be approximately at the location of the ^DLTRAP 
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TABLE II. OOT SEP Facility Test Deposition Rates (tig/pulse) 














wall when the thruster is installed in the facility-. As will be seen 

in the following sections, the measured backscatter in the MDtSINK facility, 
for otherwise identical conditions, is roughly a factor of five less than the 
^Vlar target backscatter, indicating that the cryogenic-anechoic walls of the 
MOLTRAP do provide an improvement in the backscatter levels of the PPT thruster 
discharge. 

Figure 20 shows a plot of QCM signal versus the tank target tilt angle for 
the various QQI radial locations. As can be seen, the data taken at the larger 
radii depends strongly on this tilt angle, indicating that, except for the 
data taken at 48 qn, the tank backscatter is affecting the QCM signals. This 
data also indicates that a target position of 45° minimizes the QCM 
signals and hence the wall backscatter. Figure 21 shows the QCM signals 
plotted versus radius for a target position of 45°. As can be seen, for 
smaller radii, where the plume-wall backscatter is presumably a minimum, 
the data drops with increasing radius. A comparison between the large and 
small collimator aperture angle data indicates the expected drop in signal as 
the aperture angle is reduced ; however, the small aperture ^gle data is still 
well above Of the QCM resolution, indicating that the 0° aperture can be 
used in the MDLSINK PPT thruster backflow measurements. 

3.2 MOLSINX Test Configuration The PPT thruster was mounted in the MOLSINK 
facility in a manner ide »tical to the installation of Phase I of this investigation, 
supported by a shaft entering the M3LTRAP through the upper MOLSINK doors. 

This shaft allowed the thruster to be rotated about its axis, so studies of the 
azimuthal plume behavior could be made. The thruster fired directly down to the 
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Figure 20. SEP Facility Target Backscatter 
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Figure 21. SEP Facility QCM Collimator Variation 
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lower MOLTRAP wall from the nozzle exit plane set 4.8 cm above the MOLTRAP 
horizontal midplane. The thruster power and control leads were fed into the 
^PLTRAP through the upper MOLSINK doors and down along the si^port shaft to 
the thruster. 

Eight collimated QCMs were mounted in pairs, in four rows, 38 cm, 54 cm, 

70 cm, and 86 cm from the thruster axis, respectively. A pair of QCMs Was 

used at each radius to provide some redundancy in case of failure and to increase 
the measurement accuracy, Each pair was mounted perpendicular to an cooled 

support pipe running radially outward from the thruster axis. A diagram of one 
collimated QGM pair is shown in Figure 22. The QCM collecting electrodes were 
placed 16.0 cm apart to leave room for the curved front faces of the collimators. 
A 0.8 cm square hole was cut in the QCM faceplates which were mounted rigidly 

to the QCMs, The collimators rotate about pivot points at their outer ends 

! 

and are controlled via a linkage to the outside of the MDLSINK tank. The 
aperture angle Was set at 0°; hence the total viewing angle including the 
penumbra is roughly 12°. The QCMs were rigidly mounted to the support 
pipe to prevent problems with movement of the electrical leads at low 
temperature and to provide adequate thermal conduction to the central LN2 cooled 
support shaft. Because the normal to the QCM surface was fixed, while the 
collimator axis was free to- rotate to various dip angles, a correction to the 

: I 

measured data is necessaiy. This correction consists of multiplying the mea- 

! 

sured QCM fluxes by the cosine of the angle between the QCM normal and the 
collimator axis, and accounts for the change in QCM collecting area perpendicular 
to the collimator axis. 
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SIDE VIEW 



QCM ELECTRODE 
COLLIMATOR FACE PLATE 


Figure 22. MOLSINK Collunated QCM Pair 
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Figures 23 and 24 show views of the completed array taken from the side 
and bolow; and from the side and above, respectively. The collimator 
apertures can be seen as 0.8 cm wide slits in the curved front faces of each 
collimator in Figure 23. Each QCM and its associated electronics box (mounted 

on white Teflon insulators) can be seen in Figure 24. Each CXI'I is shielrled arounH 
the sides by Kapton sheet; however, for the photograph of Figure 24, the QQ^ 
on tlie far right is unshielded so that its mounting can be observed. Also in Figure 
24, the linkages connecting the collimators together to control the dip angle 
can be seen. These linkages consist of 90° pivot arms connected with 
lengths of smaller diameter steel tubing. The main link connecting the collimators 
to outside the tank can be seen in Figure 24 extending up and out along the 
large diameter support pipe to the right of the picture. Using these linkages, 
the dip angle of the collimators in the array can be set anywhere between 0° and 
60 

Figure 25 shows a schematic of the QCM array mounted in the MOLSINK 
diaml)er. Tlie plane of the array tilts downward away from the thruster at an 
angle of about 13°. Tlus tilt is required so the outer radii collimators 
can see past the inner ones at small dip angles. The entire array is mounted 
on a slip ring tied to the central tlmister support shaft above the thruster. 

Tliis allows the array to move axially with respect to the thruster over a 
range of roughly 43 cm. The axial position of the array will be 
identified by the axial position of the QCM electrode face closest to the 
thruster, at a radius of 38 cm. Tliis QCM can be positioned anywhere from 12 cm 
upstream to 31 cm downstream of the PPT nozzle exit plane. 
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Fi<;ure 23. ^^OLSI^'K QQI ..\rray PhotograpJi 







The LN 2 cooled support pipe for the QCM array has a smaller steel tube 
installed inside it, through which LN 2 is force fed to the lower end of the 
pipe. The U ^2 returns upward out of the tank through the annular space 
between the two tubes. Throughout the following experiments, this cooling 
system was used to maintain the array QO^ temperatures at approximately -190'’C. 
During the actual data taking phase, this temperature was kept constant to 
within + 10*^C, to insure that the QCM frequency shifts were not due to temperature 
fluctuations . 

In order to accurately measure the total plume-wall backscatter over those 
areas of the NiOLSINK wall observed by the array QCMs at small dip angles, three 
additional pairs of uncollimated QCMs were mounted on brackets on the ^DLSINK 
wall. Figuro 26 shows a photograph of one such bracket, which is L- shaped with 
a V-shaped c;ut in its vertical leg for mounting on a MOLTRAP fin. The fins in the 
f40L'rPxAP run vertically from the top to the bottom of the tank; hence, wlaen installed, 
the QCMs on the bracket extend azimuthally away from the bracket's vertical leg, 
around the MDLTRAP and thruster axis. The front plate of the QCM itself has 
a 0.8 cm hole in it such that the CJCM electrode observes the volume subtended 
by a 540 half angle cone around the electrode axis. The bracket cross-section 
was designed to provide adequate conduction cooling of the QCMs to the NDLTRAP 
wall. In fact, this cooling was great enough so that the QCMs had to be heated 
with their internal temperature regulating resistors in order to maintain an 
operating temperature of -190°C. Two QCMs were installed on each bracket to pro- 
vide redundancy. As shown in Figure 25, the brackets were installed on the 
MOLTRAP wall at angles of 45°, 60° and 75° from the MOLTRAP axis, reverenced 
from the tank center point. 
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In order to insure adequate tenperature cmtrol of the test assembly, 
a number of thermocouples were placed in sensitive locations throughout the 
MOLTRAP volume. Each QCM and its associated electronics package was 
individually monitored, and if necessary, ten^erature regulated with feedback 
controllers. The internal tenq>erature of the thruster was continually control- 
led and maintained at 20-26°C at all times to prevent the oil-filled ca- 
pacitors from freezing. Finally, the lower MDLTRAP door was monitored to 
insure that tlie PPT plume did not materially affect the temperature at this 
location. 

3.3 NDLSINK Backflow Test History The PPT plune backflow measurements 
using the previously described test set-up spanned a period of approximately 
three months, with over 700 hours of accumulated facility operation. A 
typical test sequence started with sealing the outer MOLSINK doors and pump- 
ing down both the inner and outer vacuum chambers to approximately 10 ^ torr. 
During the entire process, the QCM electronics and the various thermocouples 
were continually monitored for acceptable operation. Once a low enough 
pressure was established, the LN 2 cooling of the guard vacuum walls was 
begun, and the inner chamber was isolated from the outer one. Finally, the 

glle flow was started, and the facility was allowed to come to thermal equilibrium 

-12 

at an inner chamber pressure of about 10 torr. Hie thruster was then 

started at a nominal rate of one pulse every 20 seconds, and the facility was 

allowed to equilibrate again at an average pressure somewhat greater than 
- 1 ”’ 

10 torr. Attempts at jueasuring this average pressure using a vacuum 
discharge gauge failed due to the PPT discharge interference ; however an 
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upper bound on this pressure was detennined to be about 10 torr. Once the 
entire test set-up reached this equilibrium pressure and ten^jerature, data 
was taken with the QCM diagnostics. Except for interruptions due to 
mechanical problems with the test set-up or the MDLSINK facility, the test 
was run continuously until sufficient data at all axial positions and collimator 
dip angles was accumulated. 

Prior to the first backflow measurements and with the PPT thruster 
not operating, the output frequency stability of the test array Q04s was 
measured. With the thruster not operating, the QCM mass accumulation rates 
are zero; hence the output frequency should be constant with time except 
for drift due to temperature variations. This drift was monitored for a 
six hour period fuid was found to average less than 0.6 llz for all the test 
QCMs. Tlie worst drift was found to be 2.0 Hz; hence to insure the accuracy 
of the QOI mass accumulation measurements , the total frequency shift for each 
backflow data point should be greater than 10 times this value or about 20 Hz. 
Hie collected backflow data was measured over an exposure time sufficient to 
accumulate this minimum frequency shift, except when these times became 
impractical ly long. 

Using the collimated QCM array, backflow data was taken at three axial 
locations, as measured by the axial position of the QCMs on the array closest 
to the thruster axis. These QOIs were positioned at 11.1 cm upstream, 2.54 cm 
downstream, and 30.5 cm downstream of the PPT nozzle exhaust plane. At each 
of the axial locations, data was taken at various dip angles from 0° to 
60°. Some typical QCM output frequency signals are shorn in Figures 
27 and 28 versus observation time. Both the output frequency and time are 
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referenced to zero at the begiiming of the particular obHcrvation. i'igurc 27 
shows an example of typical data taken at eitlier of the two downstream axial 
locations and for the large dip angles at the upstream location. Hie slopes 
of these data were calculated using a linear regression analysis and then 
used to calculate the QCM mass accumulation rates by multiplying them by the 
Q01 calibration constant and dividing tliem by the thruster pulse rate. The 
error in the calculated slope of the data is equal to tlio square root of one 
minus the square of the linear regression coefficient (v^l-r^). For Figure 27 
the coi'relation coefficients are around 0.999; hence the errors are very low. 
Although one of the data sets shown in Figure 28 is less accurate, the slope of 
these data is also estimated by a linear regression analysis, as before, except 
now the correlation coefficient is low and resulting error is large. 

Once the QCM mass accinnulation rates were found from tlie frequency sliift 
data, they were corrected for the difference between the collimator dip angle 
and tlie QOI surface normal, as previously discussed. Tlte final results are 
shown for the various axial positions and dip angles in Table III, along with 
their individual regression analysis correlation coefficients in parentheses. 
QfMs 1 and 2, 3 and 4, 5 and 6, and 7 and 8 are each at virtually identical 
locations and were expected to give identical results. As can be seen, tl;e 
data from the individual QCMs in these pairs can vary by as much as 50%. 
lliis variation is not consistent, but in fact, clianges with dip angle and 
axial location. Althougli tlw azimuthal separation between the QCM pairs 
is small, this may be tlie cause of the signal difference. 
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TABLE III. QCM Array Mass Fluxes (10~^ pulse"^) 


Axial Dip 
Position Angle 


30. Son 
down- 
stream 


2.54cm 

down- 

stream 


38 cm 



203. 

(.992) 

179 

(.998) 

114 

(.998) 

80.9 

(.998) 




Radius 

54 cm 


21.9 

(,984) 

38 

(.992) 

26.8 

(.996) 


70 cm 


44.0 

(.996) 


??5§8) fjsio) J:^95) 


13.8 

(.999) 


20.9 

(.999) 




8.92 

(.999) 

I'Ms) 

9.73 

(.970) 


12.5 

(.999) 

1?§^8) 

16.2 

(.982) 



.95 10.8 

(.478) (.945) 

18.4 17.3 

(.997) (.998) 

t?9i8) ^^8) 
J?5l9) ??§§9). 

I^§a6) }:§i8) 


8.0 

(.998) 

5.85 

(.998) 


7.36 

(.948) 


55.9 

(.§76) 

12.2 

(.904) 

i?4l) 

?:ao 4 ) 

?:??2) 


28.9 

9.56 

10.9 

30.9 

(.974) 

(.920) 

(.872) 

(.936) 

47.9 

29.5 

25 

35.3 

(.991) 

(.986) 

(.990) 

(.992) 

37.7 

16.7 

21.5 

31.8 

(.994) 

(.994) 

(.979) 

(.995) 

57.1 

24.1 

44.7 

43.7 

(.997) 

(.983) 

(.987) 

(.994) 

?.'Sa6) 


i’8?s) 


10.6 

3.57 

3.58 

5.63 

(.998) 

(.991) 

(.991) 

(.997) 

t-Mv 

?!9i6) 

f:§J4) 

7.5 

(.999) 

6.98 

4.58 

5.32 

9.05 

(.998) 

(.997) 

(.996) 

(.999) 


5.46 

(.995) 



9.73 

4.86 

7.36 

14.9 

(.966) 

(.894) 

(.866) 

(.986) 


18.9 

(.965) 

33.8 
(.990) 

26.9 
(.990) 

32.2 

(.973) 


5.61 

(.996) 


8.71 

(.999) 

10,4_. 



5.36 

(.997) 

10.6 

(.984) 

12.4 

(.959) 

I'Mb) 


t?9?6) 

15.7 

(.997) 

11.0 

(.998) 

12.0 

(.996) 


H§79) 


10.8 
(.960) , 

1.49 

(.051) 

5.72 ! 

(.348) I 

1 



(.§53) 

4.42 

(.718) 

?:i§5) 

tiilo) 



8.07 

(.972) 


8.87 

(.877) 


?:a§ 7 ) t? 9 fe) 


8.76 

(.991) 


8.65 

(.986) 


|!§96) I (9^4) 

14.9 ’16.7 

(.806) 



To chock this possibility md to detennine if the measured backflow has any 
significant azimuthal dependence, data was taken for two thruster azimuthal 
locations separated by 90° for the upstreajn array axial position, 
as shown in Table III. Figure 29 shows the data taken from the QCNI pair 
at the 38 an radius plotted versus dip angle for the two thrustei azimuthal 
positions. As can be seen, the different data points do not vaiy in a consistent 
manner, indicating that the scatter is primarily due to random variations in the 
measurements. Ihe data taken at the remaining three radii, 54, 70, and 86 cm, 
behave in a similar manner and indicate that within the measurement error, 
the backflow is azimuthally uniform. Hence, the variation in data between 
each Q(M of a given pair will be taken as due to random error, and the two 
signals will be averaged for the upcoming backflow analysis of Section 4.0. 

IXiring the backflow measurements of Table III, the plisne-wall backscatter 
was monitored continually with the three QCM brackets shown in Figure 25. 

The data from these QCMs was found to be independent of the arra)" axial 
position and QCM collimator dip angle, as it should be. More importantly, the 
data was also found to be independent of the tliruster azimuthal position, as 
expected from the plume s>mimetry data of Section 2.0. Tlie data taken with 
these QCMs is shown in Table IV, with respect to the bracket angular displace- 
ment from the MOLSINK center axis. 

The second QCM at the angular position of 75° failed soon after test 
inception, so only one datum is available at this location. The data 
at the other two locations shows a self-consistency similar to that found with 
the collimated array QCMs, tmd thus will simply be averaged at each location 
for the final data analysis to follow. 
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TABLE IV. Backscatter QCM Mass Fluxes 


QCM bracket angle 

1 ’ ^ Zo * 

1 Mass Flux (vig-cm -pulse ) 

45° 

' * ■ 14''~ 

j 5,18 X 10 ^ 

i 45° 

^ 5.46 X 10"'^ 

0 

0 

i 2.97 X 10'^' 

i 60° 

I 

1.51x10"^ 

1 

i 75° 

1 

1.09 X 10"'^ 


5.4 Sunimaiy Based on the experience gained during the Phase I segment 
of the PPT plume character^.zation^^^ and following a series of tests in the 
SEP vacuum facility, an array of collimated QCMs was designed and built to 
measure the PPT pliime backflow. The quality of the design was evident in its 
trouble free operation at LN 2 tanperatures and in the relative accuracy and 
consistency of the output data. Although this data behaves in a manner some- 
what different than was expected (conqjare Figure 29 to Figure 1) , the error 
bars on the data are small enough to determine approximate signal variations 
with array axial position and collimator dip angle (see Appendix 1 ). In 
general, these variations indicate that at greater axial distances upstream 
of, and at greater radial distances away from the thn:ister nozzle exit area, 
the QCM signals decrease. In addition, these signals generally increase with 
dip angle at differing rates, presumably depending on the relative dominance 
of the plume backflow or the plume-wall backscatter. Further interpretation 
of the data must await the more detailed analysis of the next section to more 
fully distinguish between these two sources of QCM signal. 
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4.0 backflow test data analysis and modeling 


During the design and operation of the collimated QCMs described in the 
previous section, the ultimate requirement for a relatively complex analytical 
reduction of the data ms always considered. Tlie use of a collimator to separate 
out the effects of the plume -wall backscatter leads to the requirement that, for 
useful results, the geometric effects of tlie collimator must be removed from 
the measured data. In addition, the contribution of the plume-wall backscatter 
must be estimated to insure adequate resolution of the actual plume backflow. 

Using this analytically corrected data, the total backflow flux through a 
representative area was calculated by integrating the data over the collimator 
dip angle. Finally, an attempt was made to reduce the data to the form 
of a scattering source function in the thruster plume. It was hoped that tliis 
source fiuiction could he used to extrapolate the calculated backflow fluxes to 
regioiis outside of the measurement area. 

4.1 Wall Backscatter Correction Since it has been concluded that scattering 
from the plume has near azimuthal symmetry, the average signals from the 
side-by-side mounted QCM pairs are used for this analysis. The pairs are 
located at distances 38, 54, 70 and 86 an from the PPT plume centerline and 
are labeled A, B, C, and f), respectively . The position of pair A relative to the 
PPT exit plane is z , wliile succeeding pairs are each offset 3.6 cm dowistream 
relative to tlieir preceding pair (sec Appendix 2). The results of the wall backscatter 
measurements are considered by calculating upper bound corrections (assuming no 
attenuation} to the array QG4 signals. Let the position on the elliptical tank 
wall be given by the angle f from the tank and PPT centerline (Figure 30). 





The wall scattering is assumed to vary as tlie cosine of the angle, 0 , from the 
\\rall normal; reduction of wall QCM signals to intensities is discussed in 

i 

Appendix 3. The wall Q(3fe were placed at y=45o,60°, and 75 ; and the measured 
results arc presented in Figure 31. The normal intensity varies exponentially with 
respect to the angle i‘. (The plume center value indicated by the uitersection of the 
straight line of Figure 31 with 4'=0 is consistent witii the Phase I measurement.) The j 

I 

wall intensity as given in the figure is used as an input source to calculate the array ; 
QCM signal contributions assuming that all the particles leaving the wall in the 
direction of a given QtM reach it. Details of the effect of QCM geometry and location 
on this calculation are fouiid in Appendix 4. The array QCM pair data and upper bound 
\\rall backscatter corrections are presented in Table V. For dip angles greater 
tlian 40°, the upper bound correction is considerably larger than the QQ^l signal. 

This means that there must be considerable attenuation of the backscattered 
wall flux (by nearly an order of magnitude) and that collisional effects in the 
plume are important. Since the attenuation is an uriknown, that part of the 
QCM signal due to plume backflow (total signal minus wall backscatter) cannot 
be known for tlie larger dip angles. Only QCM signals whose upper bound corrections 
are comparable to or less than the signal may be treated as plume backflow. 

This reduces the total number of useful data by about 20“. 

4.2 Backflow Flux Collimator Correction and Integration The QCM signal results 
from collection of particles over a solid angle defined by tJie aperture of the 
collimator and over an area defined by the opening in the back of the collimator. 

Division of signals (in mass rate) by the solid angle-area product gives their 
intensities. Total flux through a reference plane can be estimated by integration 
over angle of the product of intensity and the cosine of the angle between the 
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TABU', V. Q(’M Array Mass Flaxes and Backscatter Corrections 


Data Point 


Dip AiTglo C) 

QCM Pair 

— 
■ tJg 

QCM Signal (Pulse) 

JlE- 

Correction (Pulse) 

1 

“11. 



1 

18 

A 

7.18E-5 

6.55E-6 

2 



18 

B 

5.89E-5 

8.25E-6 




18 

C 

1.78E-5 

l.OSE-5 

4 



18 

D 

4.56E-5 

1.30E-5 

5 



26 

A 

9.27E-5 

1.15E-5 

6 



26 

B 

5.44R-5 

1.53E-5 

7 



26 

C 

2.77E-5 

1.97E-5 

8 



26 

D 

3.31E-5 

2.52E-5 

9 



36 

A 

7.54E-5 

2.25E-5 

10 



36 

B 

3.97E-5 

3.02E-5 

11 



36 

C 

2.61E-5 

3.95E-5 

12 



36 

D 

4.56E-5 

5.07E“5 

13 



60 

A 

5.89E'5 

7.76E-5 

14 



60 

B 

3.79E-5 

1.00F.-4 

15 



60 

C 

2.61E-5 

1.25E-4 

16 

T 

60 

D 

5.56E-5 

1.35E-4 

17 

2. 

54 

'0 

A 

8.91E-5 

2.06E-6 

18 



0 

B 

2.00E-5 

2.76E-6 

19 



0 

C 

6.53E-6 

3.29E-6 

20 



0 

D 

1.31E-5 

3.90E-6 

21 



12 

A 

6.85E-5 

6.62E-6 

22 



12 

B 

3.6SE-5 

7.96E-6 

23 



12 

C 

1.41E-5 

9.70E-6 




12 

D 

2.22E-5 

1.17E-5 




24 

A 

5.48E-5 

1.49E-5 




24 

B 

2.71E-5 

1.94E-5 




24 

C 

1.90R~5 

2.45E-5 




36 

A 

4.19E-5 

3.13E-5 

30 



36 

P 

2.52E-5 

4, HE- 5 

31 



36 

C 

1.9SE-5 

5.30E-5 

32 



36 

D 

3.50E-5 

6.70B-5 

33 



48 

A 

3.59E-5 

5.79E-S 

34 



48 

B 

2.76E-5 

7.60E-5 

35 



48 

C 

2.24E-5 

9.71E-S 

36 



48 

D 

4.15E-5 

1.20E-4 

37 



60 

A 

4.19E-5 

9.33E-S 

38 



60 

B 

2.76E~S 

1.19E-4 

39 



60 

C 

1.98E-5 

1.4SE-4 

40 



60 


5.04E-5 

1.42E-4 




















TABLE V. (cont'd) 


lUta Point 

r - 

ZaCciiO 

Dip Angle (°) 

QCM Pair 

QCM Signal (lAilse) 

- - 

UH 

Correction (Pulse) 

41 

m 

5 

0 

A 

7.74E-4 

6.92E-6 

42 



0 

B 

9.40E-5 

7.64E-6 

43 



0 

C 

3.79E-5 

8.58E-6 

44 



0 


9.23E-5 

9.73E-6 

45 



16 

A 

7.18E-4 

2.01E-5 

46 



16 

B 

1.72E-4 

2.45E-5 

47 



16 

C 

1.09E-4 

2.97E-5 

48 



16 

D 

1.38E-4 

3.57E-5 

49 



24 

A 

4.72E-4 

3.19E-5 

50 



24 

B 

1.30E-4 

4.01E-5 

51 



24 

C 

7.70E-5 

4.95E-5 

52 



24 

D 

1.18E-4 

6.03E-5 

53 



40 

A 

3.14E-4 

6.94E-5 

54 



40 

B 

1.95E-4 

8.87E-5 

55 



40 

C 

1.33E-4 

l.llE-4 

56 



40 

D 

1.47E-4 

1.36E-4 








intensity direction and plane normal. This procedure is detailed in the 
following discussion. Since intensity variations transverse to the plmnc are 
not determined by use of the slit collimator geometr)', only variations with 
respect to the dip angle can be calculated. Details are given in Appendix 5, 

The slit admission angle is 0.1 radian or 5.73°. Ti'ansvcrse-iiiean 
intensities (vjg-pulse -an -rad ) are given in Table VI CQCJ'f location and 
dip angle may be found in Table V). Nomiial fluxes through a plane parallel to the 
Q(H array holder which tilts at an angle of 12.7” are estbiiated (Figure 32). The 
angle a=(12.7 + dip angle) is that from the plane surface and sin a is equal 
to the cosine of the angle from the normal, thus the partial flux 

F(a) * f 1(a) sina da 

Jo 

where I is obtained from Table VI. Plots of the integrand and a similar expres- 
sion with the upper bound wall backscatter correction included are given in 
Figure 33 for a typical QCM position. For this case, the uncertainties in 
net plume backscatter do not permit integration beyond an angle of about 
450 . Partial fluxes and the associated values of a are listed in Table VII. 

Values were obtained by graphical integration. It was noted that in the range 

2 

25 ° ^ a ^ 50°, partial flux F varied approximately as for those cases of 

larger This behavior suggests that the partial integrals for small 

may be extrapolated to larger values. A least square curve fit of the form 


with axial variations removed through division by the r=54 cm values gives 
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TABLE VI. Backflow Intensities 


Data Point 

-1 -2 j-1a 

Intensity (yg-pulse -cm -rad ) 

1 

1.79E-3 

2 

1.47E-3 

3 

4.43E-4 

4 

1.13E-3 

5 

2.30E-3 

6 

1,35E"3 

7 

6.88E-4 

8 

8.22E- 

9 

1.92E-3 

10 

l.OlE-3 

11 

6.66E-4 

12 

1.16E-3 

17 

2.39E-3 

18 

5.37E-4 

19 

1.75E-4 

20 

3.520-4 

21 

1.73E-3 

22 

9.28E-4 

23 

3.56E-4 

24 

5.60E-4 

25 

1.36E-3 

26 

6.72E-4 

27 

4.71E-4 

28 

6.43E-4 

29 

1.07E-3 

30 

6.43E-4 

31 

4.97R-4 

32 

8.93E-4 

41 

2.08E-2 

42 

2.52E-3 

43 

1.02E-3 

44 

2.48E-3 

45 

1.79E-2 

46 

4.29E-3 

47 

2.72E-3 

48 

3.44E-3 

49 

1.17E-2 

50 

3.23E-3 

51 

1.91E-3 

52 

2.93E-3 

53 

8.17E-3 

54 

5.07E-3 

55 

3.46E-3 

56 

3.82E-3 
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PPT 




Figm'e 32. Backflow Integration 




i'iyuro 53 . Backflow Int 054 rand 



TABLE Vn. Plume Backflow Results (X lO'^yg-cm'^-pulse"^) 


Array 

Axial 

Position 

1 

Lad ins' 

38 cm 

54 cm 

70 an 

86 cm 

30.S cm 
Downstream 

54 

(54°) 

15 

(54°) 

3.4 

(36°) 

5.1 

(36°) 

2.5 an 

4.2 

1.2 

.12 

.33 

Downstream 

(44°) 

(34°) 

(22°) 

(26°) 

11.1 an 

6.3 

4,4 

.83 

1.4 

Upstream 

(54°) 

(45°) 

(30°) 

(36°) 


(The flux at r-38 an for the upstream position is probably low due to partial 
screening by the PCT) . Using p=2, the coefficient kCz) can be estimated. 
Results are (a in degrees) : 


Axial Position 

ci 

(cm) 

" 

k 

(pg-pulse” ^ -degree ' ^) 

St iuulard Deviat ion 

-1 -2 
(yg -pulse -degree ) 

30.5 

2.09 X 10"^ 

8.2 X lO"*^ 

2.5 

2.75 x lO"'’ 

1.0 X 10“'^ 

-11.1 

5.49 X 10'^ 

2.1 X 10“^ 


For each position the RMS error is slightly below 40%. The rough fit gives 
consistently high fluxes at r=70 an and consistently low fluxes at r= 8a. Using 
tlie above k to remove the mean axial dependence of the partial fluxes and using 
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2 

the assuned a variation to shift the fluxes to a„„ *S0° values, the radial 

max 

dependence and data scatter are shown in Figure 34. Although another radial 
dependence could give a slightly better fit, it is felt that inaccuracies in 
the data and in the extrapolation procedure do not make a more stringent curve 
fit necessary or desirable. 

Tl\e plune expansion angle is approximately 35° , thus the direction 
parallel to tlie plume boundary is approximately a»140° , giving I«0 for a > 140°. 
It may be expected that the intensity drops to small values at someviiat smaller 
values of a and that its maximum is somevdiere in the region a = 50° or slightly 
larger. Intensity plots support this conclusion. Total fluxes may be estimated 
by multiplying the 50° values by a factor of two or three. 

4.3 Source Function Studies A limited attempt has been made to numerically 
estimate the plume source function needed to duplicate the measured intensities. 

A simple model based on a linear combination of source elements is assumed. The 
QCM readings can then be represented as a known vector, equal to the product of 
an influence matrix containing QOl and source geometry effects and an unknown 
source coefficient vector. The source vector can be found by a powerful general- 
ized matrix inversion technique known as singular value analysis. Unfortunately, 
the linear source model was found to be inadequate in that there was a marked 
tendency towards partially negative sources. Also, a volume source distribution 
in the plume was, at best, not very accurate. (This is not too surprising since 
the model does not explicitly consider attenuation, which has been found to 
be important for wall scattering fluxes. 3 Nevertheless, the attempt prod\iced 
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some qualitative information and indicated the type of approach necessary for 
a possible source calculation. The model and its results are described briefly. 

A given source element (volume or surface) may be totally unscreened 
from tlie collector opening by the QCM aperture slit, or may be partially screened, 
or may be totally screened from the QQA collector. The corresponding viewing 
regions of the QCM collimator are labeled illuminate, penumbra and umbra (see 
Figure 16). For the second case of the penumbra, the collecting area is the 
overlap of the collector opening and aperture projection on the QQM plane, 

Figure 35. Calculation of the area involves considerable algebra; details 
are given in Appendix 4. The source element is assumed to contribute to a QCM 
signal an amount proportional to the element strengtli times tlie solid angle of 
the collecting area relative to the element. The first effort assumed constant, 
isotropic sources in rings of given radial and axial intervals (for example. 

Figure 36 shows a ci-oss- sect ion of 30 such rings). M influence matrix element 
describes the not effect of a given ring of imit strength on a particular QQM. 

This is calculated by dividing tlie ring into small pieces of sice Ar x AO x Az, 
calculating the collecting solid angle of the centroid of each piece, and summing the 
calculations over all volume pieces. Denote the QQ^I signal by R, source 
strength by S, solid angle by n, volume by V, then 

% ' J My S., iy. 

My = fdu. dV. = AVj. 
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Figure 36. Ring Volume Source Geometry 







The S-vector is to be found as a least square fit, i.e., to minimize the residue 

2 

defined as ?:• (R. - K- M. .S.) . IMs is done by decomposing matrijc M into 
^ ■* J J 

[M] » [EI3 [Dl [V]"^ 

T T 

where U and V are orthogonal matrices of eigenvectors of [Ml [Ml and ["Ml [Ml, 
respectively, and D is a diagonal matrix of nonincreasing quaiitities knom as 
the singular values of M. Then 

(S) = m [Pr^[Ul'^(R) 

is easily found. The ratio of the first (largest) singular value to smallest 
nonzero value is the condition number for M. If the log of this number is 
larger than the number of significant digits in the input matrix, the smallest 
value represents numerical noise and should be nulled. A study of the behavior 
of the residues of the sequence of soluticffis obtained by successively nulling 
ever larger singular values enables one to make judgments on the accuracy of 
the input and to select the proper solution vector from this sequence. Singular 
value analysis can be used to handle ill-conditioned matrices; this capability 
was found necessary for this particular problem. 

The calculated source functions were invariably negative for the larger 
radii rings in the plume. Best results were for narrow (plume angle 15° or less) 
sources with only an axial variation and with QCW signals noimalized to unity 
(relative fit) . Minimum deviations were larger than 100% RMS. This is probably 
due to attenuation effects in the actual situation. 

Since attenuation seemed important, a second effort to calculate a source 
distribution was made by assuming a surface source near the PPT plume edge 
(see Figure 37) . The source flux vector was defined in a coordinate system 
made up of a surface normal vector, n, a surface vector through the source cone 


- 91 - 



vertex, p, and a surface tangent vector, q, in the n x p direction. Quadratic 
variation with respect to axial distance and direction angles C©,i|>) was assumed. 
Although a fit with ns low as 25% l^MS error was possible, again the problem 
of a partially negative source distribution occurred. 

The results indicate that if any further attenpt at a source calculation 

l>e made, then a nonlinear procedure should be adopted, A surface source function 
should be chosen dependent on position, direction, angles, and paraaieters a, b, c, . 

such that it is nonnegative for all values of the parameters (nonlinear in 
parameters). QCM signals result by integration over solid angle and surface. 

The value of the paranieters may be found by finding a minima of the residue 
using Ne\rton's method. The difficulties that may be ejqpected with this method are: 
I'l) Selection of the form of the source function. There is little 
physical basis for a selection procedure. 

(2) Selection of initial values of the parameters. Particular initial 
values may yield a local minima but not a global one. 

(3) Stability of the calculation. 

4.4 Summary The relative con^lexity of the previous analysis is due in part 
to the plaTined sacrifice of analytical simplicity in favor of experimental 
simplicity. Several experimental design features were corrected for in this 
analysis, including the use of slotted two-dimensional collimators and fixed 
QGMs iirounted separately from the collimators. Despite these complexities and the 
error bars on the experimental data, a reasonable estimate of the total plume back- 
flow from the plume region close to the thruster was obtained. Throughout the 
analysis, conservative assumptions were made where necessary, in order to arrive 
at an estimate of tlie backflow which, at worst, is too large. Although the attempts 
at modeling this backflow in terms of a distributed source met with little success, 
a possible method was identified which may prove feasible with further study. 
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5.0 NOZZLE DESIGN STUDY 


The existing rectangular ceramic nozzle on the PPT was designed to 
help control the radiated EMI frcwn the discharge and to minimally interfere 
with the plume flow^^^\ The results of the plume studies of section 2.0 
suggest that the PPT plune has a large component of neutral species, which 
would be unaffected by the electric and magnetic fields of the discharge. 

TJie expansion of this neutral plume con^onent downstream of the discharge 
chamber may be strongly affected by the nozzle design; and hence an appropriate 
nozzle may serve to reduce the neutral plume backflow. To deterndne if this 
hypothesis is correct, a new nozzle was designed, installed on the thruster, 
and tested in the ^^LSINK facility with the collimated array for any 
changes in the plume backflow between it and the original nozzle. 

One of the primary difficulties with using QOis to measure the mass flux 
rates in the PPT plume is in their inability to resolve the plume mass flux 
versus time during a single disdiarge pulse. In addition, some questions 
exist as to the accuracy of the QCM measurement (see Appendix 1) ; and lienee 
it would be desirable to have an alternate method of measuring the plume mass 
fluxes as a way of corroborating the QCM results. One possible method which may 
prove useful is the use of Faraday cups to measure the charged particle flux 
in the plume. During the testing of the new nozzle in the ^DLSINK facility, 
a Faraday cup was installed, and its usefulness and accuracy in measuring the PPT 
plune were assessed. 

5.1 Nozzle Design, Installation and Test The original nozzle on the PPT 
expands at a half-angle of about 15” to a final exit area of approximately 
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11.5 X 16.5 cm. A straight segment of about 1 cm in effective length is 
mounted at the plume exit to provide mechanical »'inforcement of the nozzle 
joints idiich are simply epoxied together. To minimize the ablation of the 
nozzle surfaces, the nozzle was fabricated from a high temperature ceramic 
fMykroy) . 

Based on the measured plume expansion angle of 30-40 ” (soo section 
2.0), the new nozzle was designed to expand at a half -angle of 30 ^ in 
both the vertical and horizontal directions out of the PPT discharge chamber. 
This new nozzie will hereafter be referred to as the 30° nozzle. The final 
exit area of this nozzle is a rectangle 19.6 cm wide by 25.6 au high, giving 
it an area roughly twice tliat of tlie original nozzle. This larger area was 
intended to expand the plume neutral component to a pressure more nearly equal 
to the vacuum environment , and hence reduce the backflow around the nozzle 
exit. Figure 38 shows two cutaway views of the 30° nozzle drawn with solid 
lines and superimposed on the original nozzle drawn with dashed lines. The 
side view is a cutaway view in the plane containing the plume axis and a per- 
pendicular line coiinecting the electi’odes, while the top view is a cutaway 
view in the plane containing the plume axis and a perpendicular line connect- 
ing the side-fed Teflon propellant bars. To further reduce the backflow 
with the 30“ nozzle, a flat plate shield was attached to the outer lip of the 
nozzle. This plate extends outward to the dimensions of tlie thruster aluminum 
enclosure, approximately 38 cm square. Figure 39 shows the 30*^ nozzle and 
shield installed on the thruster with supports on each comer. Figure 40 shows 



a) CUTAWAY SIDE VIEW 


TEFLON 



b) CUTAWAY TOP VIEW 


Figure 38. PPT Nozzle Schematics 
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a view looking obliquely upstream into tlie nozzle throat. 

The thruster with the 30° nozzle was positioned in the MQLSINK facility 
witli the nozzle exit plane at the same location as was tlie original nozzle 
exit plane in the tests of Section 3.0. The collimated QQ4 array was positioned 
at the 2.54 cm downstreajn location and the il.l cm upstream position as in the 
original nozzle tests, and data was taken for simila* dip angles from 0 to 60°, 

Data was also taken from the mounted on tlie brackets on tlie 
MtlLTRAP walls (see Figure 25). All experimental conditions were maintained 
as close to those of the original noszle tests as possible, to provide as accurate 
a comparison becweon the two nozzles as possible. Data was taken over a two week 
period during which over 260 hours of facility operation were accumulated. 

Table VIII shows the results of the collimated QCM array measurements tabulated 
for the various array positions and dip angles. As in Table III, the correlation 
coefficients for the data are shoivn in parentheses and provide a measure of 
the error in each datum. 

The results of the plume -wall backscatter measurement from the QCMs 
mounted on the ^DLS1NK wall are tabulated in Table IX. Following an analysis 
identical to that in Section 4.1 for the original nozzle backscatter, these 
backscatter mass fluxes were reduced to intensities and are shown plotted versus 
angular location on tlie ^^0LSINK wall in Figure 41. Also shown is the data from 
Figure 31 for the original nozzle. Altliough the two sets of data are within a 
factor of two of each othei'^ and hence, are within the QC>1 error bar, the 
consistently higher data for the 30° nozzle suggests a greater mass flux 
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arrival rate at the walls at higher angular locations. This, in turn, suggests 
that with the nozzle, the average plume expansion angle is larger than 
with the original nozzle, or alternatively that the 30° nozzle -shield 
coinbitiation is directing more of the plume material downstream, and thus, de- 
creasing the plume backflow. Distinguishing between these possibilities would 
require further testing; however, some indications that the sec» nd alternative 
is correct will be seen in the following section. 

TABLli VI II. 30^' Nozzle Array Mass Fluxes (X 10 ^vig"Cm ^-pulse 




Radius 

Axial 

■ffll 

38 cm 

54 

cm 

70 

cm 

86 cm 

Position 



-2 



"4 ■ 

5 

6 

7 

8 

2,54cm 

dowisti’oam 

0° 

18.5 

f.999) 

23.7 

(.997) 

■HM 

9 2 
G946) 

11.4 

(.930) 


7.0 

(.979) 

6.5 

(.938) 


12° 

19.8 

(.999) 

20.4 

(.999) 


12.7 

(.991) 

7.2 

(.986) 


7.1 

(.954) 



24° 

19.6 

(.999) 

17.3 

(.999) 

7,0 

(.947) 

9.1 

(.915) 

6.3 

(.993) 

- 

9.1 

(.993) 

10.6 

(.986) 

» 

36° 

15.9 

(.998) 

12.5 

(.997) 

- 

13.3 

(.973) 

5.8 

(.969) 

- 

15.7 

(.981) 

17.8 

(.958) 


60° 

10,4 

(.997) 

8.7 

(.997) 


15.1 

(.967) 

9.3 

(.985) 


13.3 

(.983) 

12.3 

(.968) 

11. Ian 
upstream 

12° 

9.6 

(.998) 

6.5 

(.995) 

8.9 

(.767) 

9.9 

(.972) 

12.0 
(.984) ■ 

8.6 

(.848) 

3.8 

(.794) 

- 


26° 

14.6 

(.998) 

10.0 

(.998) 

- 

9.6 

(.925) 

6.5 

(.993) 

6.0 

(.825) 

8.9 

(.961) 

10.4 

(.948) 


36° 

9 4 
(!999) 

6.7 

(.996) 

7.5 

(.908) 

8.7 

(.861) 

7.3 

(.990) 

5.9 

(.949) 

9.5 

(.981) 

9.9 

(.962) 


60° 

10,0 

(.996) 

6.5 

(.994) 

■ 

13.3 

(.965) 

8.9 

(.989) 

4.6 

(.602) 

13.9 

(.950) 

13.5 

(.950) 
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TARI.U IX. 30'^ Nozzle Backscatter Mass Fluxes 


QCM Bracket Angle 



Mass Flux (vg-on -pulse ') 

45° 

8.49 X 10"^ 

45° 

7.13 X 10”^ 

60° 

3.88 X lO’"^ 

75° 

1.37 X 10’^ 

75° 

2.06 X 10'^ 


5.2 30^' Nozzle Backflow Analysis Using the best exponential fit of the 

measured plume -wall backscatter, as shown in Figure 41, corrections to the 
collimated QCM array data of Table VIII were calculated in a manner identical 
to that of Section 4.1. Using this corrected data to determine the maximum 
allowable dip angle, these data were integrated over the dip angle to determine 
the total backflow flux through the QCM array plane Csee Section 4.2 and Figure 
33) from the plume region between the dip angles of zero and the maximum value. 
To insure a ccmslstent comparison between the original and the 30° nozzle, 
the maxhnum dip angle for tlie individual collimated QCMs for the 30° nozzle 
was taken to be identical to that used in the original nozzle anax sis (see 
Table VJII) . Table X shows the integrated backflow fluxes versus array axial 
location and QCM radius. For comparison the backflow fluxes measured with the 
original nozzle are also slioMi along with the values of the maximum dip angle 
in parentheses. 

Thase tabulated values are plotted versus radius in Figures 42 and 43, for 
both nozzles and for the two array locations. Tlie data taken for the 2.54 cm 


- 101 - 








TABLli X. Pluinc Backflow Results for Both Nozzles (X lO'^^yg-cm'^-pulse'^) 


Array ' 
Axial 

Nozzle 

QQvi Radius . j 

Position 

38 cm 

54 cm 

70 cm 

86 cm 

2.54an 

downstream 

30° 

Nozzle 

4.9 

1.67 

0.72 

0.58 


Original 

Nozzle 

■’t‘»14 

(44°) 

1.15 

(34®) 


0.37 

(26°) 

11. Ian 
upstream 

30” 

Nozzle 

3.1 

2.2 

1.4 

0.7 


Original 

Nozzle 

6.55 

(54°) 

4.4 

(45°) 

0.9 

(30°) 

1.3 

(36“) 


downstream array axial location is, within the QCJi error, identical for both 
nozzles. Conversely, for the 11.1 cm upstream axial location, the 30° nozzle 
has a significantly lower backflow for the lower radius locations. These results 
indicate that with the original nozzle, a significant backflow arises between 
the downstream edge of the nozzle and the radius of the closest QCM (38 cm) . 

With the 30” nozzle and shield, this portion of tlie backflow is considerably 
reduced, leading to the drop in the measured data at the 11.1 cm upstream 
position. The equivalent data at the 2.54 cm downstream location indicates 
that this deflected portion of the original nozzle backflow is directed down- 
stream and not just radially outrard to the edge of the shield and back up- 
stream. Furthermore, the equivalent data at tJie 2.54 cm downstream location 
suggests that the 30 ” nozzle does not materially change the downstream plume 
flow profiles, and indicates that the shield may have the greatest effect on the 
backflow. 




















BACKFLOW FLUX, x 10 ^ fj.g~cm -pulse 



f'ii’urc 42, Rackflow Radial Variation - 2.54 an Doivnstreani 
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S.3 Faraday Cup Assessment A Farada)' cup collects and moasurcs the diarged 
particle flux to a collecting electrode from tlie jUasnia imder study. This 
collecting electrode sits at the bottom of an insulating cup and is connected 
to a large reference electrode directly in touch with the plasma, A screen 
grid, placed over the opening of tlie cup, is normally biased negatively with 
respect to the reference electrode to repel the electron flu>' to the collecting 
electrode. The ion flux to the collecting electrode flows through the circuit 
and buck to tlic plasma via the reference electrode. The ion flux or current 
is measured in the circuit tuid used to estimate the total mass flux via the 
ionization fraction. The time response of the l*arada>’ cup is limited by tlie 
electronics used to measure the ion current, and hence can be made fast enough 
to allow resolution of the mass flux during just one IW discharge pulse. 

The j)urj)Oso of this nivest igat ion was to assess (he ability of 

tile Faraday cup to provide accurate measurements of the IW plimie backflow, 

and to corroborate the previous QCM measurements. To this end, a Faraday cup 

2 

with a collecting area of 1.27 cm was mounted on the ^DLSINK tank wall in the 
PPT nozzle exit plane on the opposite side of the tlirustcr from the collimated 
QClvi array. Tlie collecting cup was turned to face directly towards the PPT 
axis, so MS to (observe a backflow mass flux similar to that measured with the 
collimated QCMs. Because it was attached to the MDLSINK wall, the Faraday cup 
operated at a temperature of about 15-20 K. 

The ion current through the cup circuit was measured using one 
of two alternate methods. In the first method, a lOOOfi resistor was placed 
in scries witli the circuit, and the voltage across this resistor was monitored. 
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Tlie second inethod used a Tectronlcs Model 0042 currant j robe which monitors 
tlie ion current through thp nii^gnetic field indticed in the circuit wiring, and 
thus, interferes minimally with the actual probe circuit. The ion current 
was monitored and recorded on gn oscilloscope triggered by the noise generated 
by the thruster discharge. Figure 44 shows two oscillographs of the ion 
current taken under identical conditions with the two alternate current measuring 
methods. The ion current measured with the resistor-voltage probe combination 
is the observed voltage divided by lOOOo , which is consistent with tlie peak 
current of about 0.8 to 1,0 ni/V, taken using tlie Tectronix curi'ent probe. As 
can be seen, tlie Tectronix probe suffers from a relatively small signal to 
noise ratio. 

In order to accurately moasuie the total ion flux to tlie Faraday 
cup, the grid bias must be set to repel the maximum number of electrons, and 
yet, minimally affect the incoming ion trajectories. To determine the 
appropriate value for the bias, measurements ivere made over a range of bias 
voltages from 0 to 100 volts, 100 volts, the gap between the grid and the cup 
would break down occasionally, thus effectively limiting the maximum bias to 
this value, Tlie mensiucd ion currents were Cciind to generally increase with 
Increasing bias voltage, ^vliile the total ion flux to the Faraday cup integrated 
over time varied by no more than a factor of two over tlie total range of 
bias values. Tliis beliavior indicates that the higher bias voltages are more desirable; 
however, within the overall error bar, the value of Mas voltage has only a 
small effect on tlie results. 
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rigure 44 »show.s that the ion current to the cup is negative during tlic 
first 30 psec following the PPT discharge initiation. Since the ciq> is approxi- 
mately 1 meter fr(xn the PPT thruster, this 30 psec delay is consistent with the 
time required for ions moving at 30,000 m/sec to reach the collecting surface 
(see Section 2.0). Following this delay, the current shows a distinct double 
peak structure which cannot be due to plasma-wall backscatter since tlie time 
separation between the two peaks is too short. A possible e)q>lanation may lie 
in the observation by Palumbo and Begun^^^^ of plasmoids 0>lobs of luminous 
plasma) being accelerated off the thruster electrodes. The separation of these 
plasmoids was roughly correlated with each reversal of the PPT discharge 
current. For the thruster mder study, the discharge current reverses once about 
15 psec from the beginning of the discharge implying one plasmoid followed 
by a continuous plasma flow. The separation between the peaks of the Faraday 
cup data is about 20 to 30 psec and may be due to the differing velocities 
of these two plasmoids driven by the decaying current. 

The tota] Ion flux collected by the Faraday cup during one PPT 

discharge pulse can be calculated by integrating the measured ion current 

over the total observation time. The average flux was obtained by estimating 

this integral for several PPT discliarges and averaging tlie results. This flux 

11 -2 -1 

was found to be 2 x 10 ions -cm -pulse , assuming no double ions. Using the 
average atomic weight of Teflon of 16.7 amu for the ion mass, the total mass 
flux collected by the Faraday cup is 5.5 x 10~” pg-cm” -pulse"^. The measured 
Q(3I mass flux corrected for the larger aperture of the Faraday cup and the radius 
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of 1 meter from the PPT thruster axis is over 20 times larger than this value, 
indicating that either the Faraday ct^ is not collecting all the plasma ions or 
that the actual plasma ionization fraction is so small that the ion mass flux 
is only a fraction of the total. In either case it appears that the Faraday 
cup does not provide an accurate quantitative measure of the actual PPT backflow 
flux. 


5.4 Summary The backflow measurements taken with the 30" nozzle on the PPT 
confirm the ability of the collimated QCM array to distinguish between changing 
mass flux levels and indicate that the overall backflow with the 30° nozzle- 
shield ccmbination is less than that of the original nozzle. It is believed 
that the major factor in reducing the backflow is the shield and not the 
increased expansion angle of the 30° nozzle. The larger plume-wall backscatter 
found with this nozzle in^lies that the mass which previously was backflow is 
being redirected to the outer edges of the primary plume, however, it nay also 
be due to a larger plume expansion angle. 


liven with the large error on the QCll backflow fliix measurements, 
the QOf still appears to be better thrn a Faraday cup for measuring this flux. 

The Faraday cup is useful in observing the time history of the PPT discharge 
plume; however, even here there remain some problems in interpreting the data. 
During these studies it was not possible to satisfactorily explain the apparent 
negative ion current during tht^ first 30 ysec of the Faraday cup signal. In 
addition, the behavior of the ion current dependence on bias voltage is dif- 
ficult to fully explain. Before the Faraday cup can be used with confidence, 
tliese points should be addressed. 
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6.0 SlIM-iAm' AND CONCUISfON 


The simplicity and reliability of micropoiuid PPT thrusters stimulated the 
development of a lay.qer millipoLuid version suitable for an expanded range 
of applications including north- south stationheeping. Tlie larger impulse 
bit and greater total impulse of this millipound thruster have led to concerns 
about potential exhaust plume contamination of sensitive spacecraft surfaces , 

.mu thus, to a program at the Jet Propulsion haboratoiy aimed at characteri;:ing 
the flow field of tlie millipound thruster plume, especially in the region upstream 
of the thruster nozzle. Phase I of this prognmi, detailed in an earlier report, 
was aimed at developing an understanding of the PPT plume-wall backscattcr 
levels in the special MOLSINK vacuimi facility, and a method of measuring the 
PPT plume backflow in the presence of this backscattered flux. 

Phase II of the PPT plume characterization is detailed in this report. 

The priman' puiiiosc of this Pliasc II effort was to measure, using the method 
devclopcil in Phase 1, the PPT plume backflow flux over a range of loc.ations 
radially away from and upstream of tiie tliruster nozzle exit plane. As a corollary 
to this effort , a secondaiy task was to develop a better understanding of the 
tliruster primary plume nuiss flux doioistroam of the nozzle exit plane. Once 
these original tasks t\fere completed, .a .small additional study was made to determine 
the .sensitivity of the PPT plume backflow to a different thruster nozzle design, 
and to the presence of a shield surrounding the nozzle. 

In order to insure the succes.sful completion of the MiLSINK facility PFf 
plume backflow iiiea.su rement, the studies of the primaiy PPT plume were carried out 
first, so that the results could Ik' factored into tiie design of the backflow 
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moil.surcmcnt system. Scvei'al tests wore made using various diagnostics 
including a ^fylar sheet deposition target, a double QQ] probe, high speed 
photography, spectroscopy and carbon disk deposition analysis. These tests 
confirmed that despite the high aspect ratio rectangular PPT discharge chamber, 
the primary plume is essentially axisymmetric at distances of over 7S cm from the 
exhaust nozzle. The pliune was consistently found to iiave an exi^ansion half 
angle of between 50" and 40", which, from the double QQI probe results, 
represents the plume radius at which the downstream mass flux has dropped 
to 10% of its centeriine value. Furthermore, the plume mass flux profile was 
found to have a half-angle at half its centerline value of about 20° » 
indicating that the plume is more collimated than was previously supposed. 
Observations of the plume with high-speed photography indicate that it 
consists of a high velocity ki/sec) luminous plasma combined with a much 
slower nonluminous gas. Mien these two mass flows are averaged, it is known 
that the average plume exhaust velocity is about 17 bi/sec, hence the slow 
nonluminous component must be a significant fraction of tlie total plume. The 
flow of this component may be less dominated by the electromagnetic forces 
which tend to confine the pluiie to the thruster axis, and thus the plume back- 
flovN? may originate with the more ordinary gasdynamic expansion of this cooler 
phmie component. 

Spectroscopic evidence confirms that the luminous component of the plume 
contains ionized fluorine and carbon, and iicnce, is highly energetic. The 
observed erosion of surfaces placed in the primary plume is most probably due 
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to impingement by this high energy plasma. Analysis of the deposits on 
carbon sampling disks indicates that the Teflon propellant plume is primarily 
vapor with few large particles. This vapor condenses uniformly over the 
collecting surface, indicating that the QCM calibration constant analysis 
is correct in assuming no point masses on the sensing electrode. The observed 
erosion of the PPT electrodes is confirmed by the presence of copper on the 
sampling disk placed in the primary plume. This copper was not observed on the 
sampling disk upstream of and just outside the nozzle, indicating that the eroded 
copper is not in the plume backflow. 

During the Phase I effort of this PPT thruster plume characterization, a 
concept of measuring the plume backflow in the MOLSINK facility was developed. 
This concept uses collimated QCMs to observe the plume backflow while avoiding 
any observations of the plume-wall backscatter. During the current Phase II, 
the design of these collimated QQis was refined based on testing in the SEP 
facility, and an array of collimated QCMs was assembled and installed in the 
MOLSINK facility. Included in this installation was a series of QQis mounted 
on the MOLSINK inner wall to observe the plume-wall backscatter. These 
observations were used to check the collimated QCM array measurements to insure 
that the backflow measurements i\rere not compromised by contributions from 
this backscatter. ilie results of these measurements were analyzed extensively 
to account for various complicating design features of the experimental system 
and to determine the actual plume backflow levels at various radial and axial 
positions in the PPT nozzle vicinity. The results indicate that in a region 
between 30 and 86 an from the thruster axis and from 11.1 cm upstream to 30.5 cm 
downstream of the nozzle exit plane, the plume backflow mass flux is of order 
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10 g -cm -pulse . Although the error bar on this flux is roughly a factor 
of two, an analysis of its radial dependence indicates an approximate inverse 
quadratic drop-off with distance from the thruster axis. 

In order to place the magnitude of the measured PPT plume backflow in context 
with other types of thruster systems, a comparison can be made to the backflow 
flux from tlie plume of an 8 cm mercury ion thruster, lliis thruster has a 
nominal thrust level of one millipound as does the PPT; however, its specific 
impulse is api^roximately 60% higher and its mass flow rate is 40% lo\\?er than the 
PPT thruster. The 8 cm mercun/- ion thruster was developed for applications similar 
to those of the pulsed plasma thruster, including stationkeeping and attitude 
control. The total backflow from the plume of this ion thruster was estimated 
by summing the contributions from both tlie mercury propellant and molybdenum 
sputtered from the grids , for locations identical to those where the PPT 
backflow measurements were made. The backflov\r flux from the PPT, as shown in 
Table VII, was corrected for the nominal pulse rate of 0.2 pps and was foimd 
to be identical, within cx|ierimental error, to that found for the ion thruster. 

'i'lic investigation of tlie primary plume indicated that a significant fraction 
of the pliune mass is relatively slow moving and at a low temperature. This 
suggests that the flow of this material is predominantly gasdynamic, and 
thus, dependent on the discliarge nozzle design. Furthermore, since the higher 
energy plasma component of the plume is confined to the magnetic nozzle set up 
by the discharge, a large fraction of the measured backflow may arise from the 
gasdynamic flow around the nozzle lip of the lower energ)'- fraction of the plume. 
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To test this possiblity, a new nozzle was designed for the PPT thruster 
with an exi^ansion angle of 30° and an exliaust plane area of roughly 
twice that of the original nozzle. Included in this design was a flat plate 
shield which extends outward, in the nozzle exit plane, to a size correspond- 
ing to the size of the thruster enclosure. 

This new nozzle-shield combination was tested in the MOLSINK facility 
under conditions identical to those of the original nozzle backflow test. 

Tlie results were analyzed and tlien compared to tliose of the original nozzle 
backflow levels. This comparison indicates that the nozzle itself makes 
little difference in the magnitude of tlie backflow, but the shield has a 
strong effect on reducing the backflow in regions upstream of the nozzle 
exit plane. This reduction is as much as a factor of three for close radii, 
while at larger radii, the backflow is not affected. 

The major pun^ose of the plume characterization was to measure the 
backflow flux from the PPT pliunc in order to assist in deteniiining the effect 
of its deposition on various spacecraft surfaces. The final assessment 
of the plimie backflow must include considerations, not only of the total 
backflow flux, but of what fraction of this flux actually deposits and what 
actual effect does tliis deposit have on the various surfaces. These problems 
depend strongly on the t>T)e of surface (solar array, thennal radiator, space- 
craft housing, observation ports, etc) and the surface temperature. For example , 

certain solar cell array designs use Teflon covers , and hence may be totally 

(1 S') 

luiaffocted by the plume deposits . At the other extreme sensitive optical 
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sensors with heavy plume deposits would have perturbed transmission or 
absorption characteristics, and thus, significantly degraded performance 
levels. Finally, continuing improvemei ,ts in the thruster performance, includ- 
ing increasing specific impulse and thrust, imply a reduction of the backflow 
mass flux. With all of these factors in mind , the current measurements can only 
be used for rough estimates of what the ultimate effect of the PPT plume 
backflow will be. 
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APPENDIX 1 


QCM ERROR ANALYSIS 

As discussed in section 2.3, the priiiiaiy PFF pliune has a significant 
neutral component. From section 5.3, it is evident that the plume backflow 
has an even larger component (up to 95%) , hence to study the PPT plume mass 
flow, a diagnostic system is needed wliich is sensitive to both charged and 
neutral particles. Quartz crystal microbalances (QCMs) were chosen for this 
study because of tills requirement and because they provide in-situ measure- 
ments which do not require removal for analysis and subsequent danger of 
contajiii nation. In addition, their relative compactness allows several to be 
used simultaneously witliout ov^ercrowding the vacuimi facility. 

The accuracy of a particular QQ^i measurement depends on several factors 
which can be split into two broad areas. The first area includes factors 
which describe the relation between the mass flux at a certain location 
and the actual collected mass when a QQi is at this location. These factors 
include considerations of the particle optics to the sensing crystal 
(i.e., collimator design and leakage), spurious mass accumulation (due to the 
vacuimi tank env i roniiicntal pressure and the pulsed tliruster operation) , and 
the accoiranodation coefficient of the collecting surface. The second area 
includes factors wliicli relate the accumulated mass to the measured output 
frci[uenc>' shift. These factors include the QCM temperature sensitivity 
and electronic stability, the value of the calibration constant, and the 
calculation of the frequency shift versus time. These factors will be 
iliscussed individually in the following paragraphs and then used to calculate 
a total error estinuitc. 
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PM^TICLE OPTICS The pluine ion number density has been measured and 

9 11-3 

estimated to be between 10-10 cm in the nozzle exit plane for distances 

of order 50 cm away from the thruster axis. Presuming that the total 

heavy particle number density is about 10 times the ion density gives 
10 12 -3 

10 -10 cm . Average heavy particle elastic collision cross sections are 
of order 10 ^^-10 cm^; hence, the particle mean free paths are around 

4 

10-10 an. For the collimator designs used in the plume study effort, 
the characteristic size is about 8 cm, which is less than the particle 
mean free path. Thus, it is consistent to use simple optics theory with 
its concepts of th^ illuminato, penianbra and lunbra viewing regions for ’■he 
collimator design and operation. Collisional effects which would smear the 
distinction between these viewing regions can be and are neglected in the 
data reduction of section 4.0. 

The possibility of mass flux leaks around the collimator to the QOl 
crystal was guarded against by careful shielding of the QCMs with Kapton 
slieet around the back and sides of the QCM body. The shields combined with 
the long particle mean free paths effectively prevented any flux from hnping- 
ing on the collecting surface except for that which entered via the collhnator 
aperture . 

SPURIOUS MASS ACCUMULATION The ultimate ^OJSINK facility pressure 
-12 

is of order 10 torr, while the temperature is of order 25 K. Assuming 
the gas in the chamber has a molecular weight equal to the average molecular 
weight of Teflon (16.6 amu), and using the ideal gas law, the particle density 
can !)c estimated as approximately 10 g-an . From kinetic theory the average mass 
flux through a unit area in one direction in a volume of gas is the product of 
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tile gas density and the sound speed divided by 4. Using 5/3 as the ratio of 
specific heats, the sound speed at 25 K is 1.44x10^ cni-sec^; thus the average 
mass flux is 3x10 g“dn^"-sec^ , This flux is far smaller than tlie measured 
QOi fluxes; hence the error introduced by this low pressure environmental 
flax is negligible. 

During the >;clual test, with the thruster firing, the MOLSINK environ- 
mental pressure is not constant, but fluctuates up to values considerably 

-12 

larger than the ultimate low pressure of 10 torr. This fluctuation is 
due to the discharge pulse mass input into the tank, and thus, is of very 
short duration. This sliort duration, combined with the pulse rate of once 
every 20 seconds, prevents the available pressure measurement Systran from 
giving an accurate value of the pressure versus time, primarily because the 
gauge time constant is about 3 seconds, i.e., much longer than the discliarge 
pulse time. Due to this lack of available experimental data, the effect 
of this pressure rise on the QOI measurement must be estimated analytically. 

The PFf discharge mass is about 1.56 mg-pulse^. This mass leaves the 
thruster and travels to tlie ^DLSINK wall, 1.5 m away, witli an average sonic 
velocity of c. iXjring tlie time after this first wall interaction, the mass 
accumulated va the QCM surfaces represents an error on tlie desired signal. 

J-rom the Phase I study, 941 of the discharge mass is absorbed by the wall 
while the remaining 6% is backscattered . The backscattered mass travels 
upstream to the opposite end of the MOLSINK, conservatively about 3 m away. 

The time required to do this is 3 meters divided by the sound speed: 3/c. 
Assuming tlint during this time, the mass density, p , is equal to tlie remaining 
(>» of tlie pliune mass divided by the tank voliuiie (6.28 iif), then the net mass 
accimiulaled on the l•eprescntat ive QCM surface is; 
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nils is the accumulated mass during the time between the first and second 
wall interactions. Assuming that all subsequent wall interactions reduce 
the net mass by 50%, then the total mass accumulated on the QO^l after one 
discharge pulse is the sum of Mj _2 ^^2-3 ^ '^^3-4 *■'* found to be 

_7 .1 

l,8xl0‘ Mg-cJir “ pulse . Tills figure is independent of the mass velocity 
and tempei'Qture; however it depends critically on the assumption of uniform 
density over the tank volume between collisions. If this assumption is 
relaxed in favor of a more concentrated plume volume, then the total mass 
accumulated on the Q(31 would drop since the time of exposure of the QGM 
surface to this density drops. In addition, if more than 50% of the mass is 
absorbed during each wall interaction, the accumulated mass would drop. 

These assumptions could easily make an order of magnitude difference in 
the analysis. 

The calculated environmental flux represents that flux which would be 

measured by an uncollimated QCM in the MOLSINK tank. Comparing this result 

with experimental data from collimated QCMs must be done after the collimated 

QCM data has been reduced to a total flux by having the collimator effect 

integrated out. This integration was done for the original PPT thruster and 

nozzle in section 4.2 and the overall results are shovm in Table VII. As can 

be seen, the smallest integrated signal, which must be greater than the 

environmental flux since it includes the backflow flux as well, is about 
-5 -2 -1 

10* yg-cm -pulse at a location of 2.5 cm doivnstream and 70 cm radius. 

Thus, the previous analytical result is roughly 180 times too large. As caii 
be seen in Table VII, the majority of the data points are significantly 
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larger than l9 ug-cm ^-pulse ^ thus the error due to the environmental flux 
is usually small. At its largest this error is no more than a factor of two 
on the smallest backflow flux value, 

QCM ACCON^IQDATION COEFFICIENT The boiling points of the basic monomers 
of depolymerized Teflon are discussed in Reference 16. The lowest boiling 
point mentioned is for Hexafluoropropene which is 144 K. Further- 

more, Reference 17 indicates that most chemical recombinations in the PPT 
plume occur within a few inches of the nozzle hence, the plume constituents 
are primarily tliese basic monomers. Since the QG4s of the current study 
were maintained at approximately 75 K, the accommodation coefficient for these 
materials should be very close to 1.0; hence all the incoming mass to 
the QQ4 surface would be condensed on it. In an}^ event, even if the accommoda- 
tion coefficient is less than one, the QCM temperature is considerably lower 
tlian ordinary spacecraft surface temperatures. This indicates that 
any spacecraft surf.acc would iiave a lower accommodation coefficient tlian the 
QCMs used here, and thus would collect less material than that collected on 
the QCMs of this study. In this respect, the measured QCM mass fluxes are 
conservative values of vdiat would actually be affecting a spacecraft. 

(jC?l THNlPFRAlimF. AND IILECTRONIC STABILITY During the Phase T segment 
of tlie phune characterization study, a design was developed wliich 
provided maxiimm' electronic stability aiui miniinimi temperature variations 
in the output freciuency. The use of a doublet temperature compensated design 
and AT cut quartz cr)*stals provides a Qll-I which is virtually Insensitive to 
temperature fluctuations in the range from -220 to -180 °C. 
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The thnister was floated with respect to the entire QCM system and the 
electronics siqjport package for each QCM was modified to use bipolar 
transistors in order to eliminate the CJCM sensitivity to thruster E^^I. 

As discussed in section 3.3, the QCM stability was experimentally checked 
prior to the plume backflow measurement. The results demonstrated that 
the maximum drift A output frequency can be up to 2 . 9 Hz . To minimize 
the error induced by this drift, where practical, each experimental datum 
was determined by at least a 20 Hz shift, thus the maximum error is less 
than 10% for most of the data. 

FREQUENCY, SHIfT CALCULATION Each Q(31 output frequency was monitored 
versus time for test times of order 20-40 hours. The time was measured 
with a clock accurate to several minutes, and the frequency shift was 
measured with a counter accurate to 1 Hz. Measurements were made on the 
average of once per hour; thus a minimum of 20 data points were accumulated 
per QCM flux measurement. A linear regression analysis was used to cal- 
culate the slope of the frequency shift versus time, which is directly 
proportional to the accumulated mass flux. The error in the slope calcula- 
tion is equal to + /T? , where a is the linear regression correlation co- 
efficient. These coefficients are tabulated with the backflow measurements 
in Tables IV and VIII. As can be seen, the usual correlation coefficient 
is about 0.95; thus the slope error is about + 10%. 

QCM CALIBRATION CONSiANT During the QCM development of the Phase I 
segment of this effort, the calibration constant was derived analytically and 
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measured e.xpcrimenfaUy. The calibration constant for an identical Q(’>1 
was ineaaurecl very accurately by using sputtered aluminum as 

the deposit material. The calibration constant depends on the quartz ciy.stal 
oscillation frequency, the cut angle, the quartz density and the uniforaity 
of the deposit on the sensing surface. These quantities are all fixed 
once the particular (^CIM design is chosen except for the deposit unifonnlty. 
This imirorinity ms examined and is discussed in section 2.4, where it is 
shenm that ti\e PPT plume deposit is veiy smooth with relatively few point 
ma.sses. Baseil upon this exiuninat ion, it was concluded that the QOIs in use 
I’or this plume study effort are Identical to that calibrated by Phillips. 

The final value of the calibration constant lised in the Phase II effort is 
the one found by Phillips, which is 1.77xl0“^ug~air‘^~ilz’’^ ± 5%. 

SUMMARY The total error of tlie QCl-l flux rnensurement is Just tiie root 
mean square sum of the errors discussed in tlic previous paragi'aphs. This 
error is im average about 1 however, for a few selected data it can be 
as higli as a factor of 2.5, due to tlie potential environmental fliLX. Once 
this data is used in tlie analysis of section 4,0, an additional uncertainty of 
about a factor of 2 is introduced due to the analytical assumptions and limita- 
tions. The final values of plume backflow mass flux Jiave on uncertainty 
about a factor of 2, up to 5 for tlie wxu'St case. 
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APPENDIX 2 


COORDINATE TRANSFORMATIONS 


A coordinate system R=(x,y,z) is attached to the thruster (see Figure 2-1). 

Tlio axial distance do\m the pliuiie from the nozzle exit plane is z; x and y are 

transverse coordinates witli the Q01 array holder located in the negative 

y-direction (x=0). The array holder slajits down and away from the thruster 

exit at an angle of 12.7'’. Four pairs of QCMs, labeled A, B, C, D are located 

on the holder at negative y-values of 38 cm, 54 cm, 70 cm and 86 an, respectively. 

Since the y increment between pairs is 16an, the corresponding z increment is 

16 tan 12,7°=3.606 cm. Each pair has a QCM positioned at x»+ 7.9375 cm. 

A translation of it to coordinates it’ = (x' y’ z') with origin at pair centers 

s s s s 

(sec pair B in Figure 2-1) is accomplished as “ it - , where 

\ = 0; y^. = -38,-54,-70,-86; and z^ = z^,z_^ + 3.606, z^ + 7.212, z^ + 10.818 for 
pairs A, B, C, D, respectively. Pair A QCMs are labeled 1 and 2\ pair B, 3 and 4. 
pair C, 5 and 6; pair I), 7 and 8. (Idd niuuber QOls have a positive x 
po.sition, even negative. 

A second transfoniiation of coordinates from it' to it » Cx^,y„,z_) consists of 

s s s s s 

an .’(.-translation Ax = 7.9375 to a particular QCM center plus a rotation about 
the x-Jixis of an angle a (I'ight hand rotation, see QCM 2 in Figure 2-1) ; 

X = X* + 

s s 

-p = sina + COSrt 

The s • axis is mnv’ (lie ()CM a>lMmator look direction for a dip angle n. 

Tlie normal to the (JCM crystal has a fixed ilip angle of 22. S"; hence the 
collector angle relative to is i = a-22,5'\ For -8‘V a then 
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Figure 2-1. Coordinates Relative to QCMs and Plume 



-30,5*^ T 34.5°. 


! 

i 

I 


If the plume is to be‘ asSiined axially symmetric, then more convenient 
coordinates are t. - (r,0,z) where x ® r cose, y - rsine. The transformation 
(r,0 ,z)»^(Xg, yg, Zg) is a three step process. 
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APPENDIX 3 


SOLID ANGLE CALCULATION FOR IVALL QCNI 


The solid aiigle of the QO-te used to measure backscatter from the wall 
is calculated neglecting wall curvature. The QCMs see the wall through a hole 
of radius R^; their surface is a distance d from the wall and h from the cover 
plate front (see Figure 3-1). The equation of the sight cone of collector surface 
element dA^ located at r * R,e = 2 = 0 is 
[r COSO - R(l- 2 /h)]^ + r^ sin^ 0 = 


Tl\e surface area seen by elaiient dA is boLind by the circle 

[r cosG - R(l-d/h)]“ + r^ sin^O = (R^ d/h)^ = R? 

\diich is centered at x = rcose Coordinates relative to the circle 

center are denoted by a prime. Consider a wall element dA located at (r,Q), 

the angle y between a line fron^ dA^ to d.A and the wall normal is given by 

d 

V(r ' COSO • -Rd/h) “ (r ' sinG ' ) ^'+d^ 


It is assLuned that the wall backscatter intensity has the form 1=1© cosy- 
Also, the solid aragle of clA^ as seen by the wall element is 
cosy dA^ 


do = 


(r ' cos 0 ' - Rd/h)“ + (r'sino')^ + d^' 


Tlie total flux to the QCM surface is 

2n R, 


>-/ 


hhi cl\ = 


I dYdA f r’dr'do* 

° J 0 Jo [(r ’ COSO' -R d/h) (r ' sine ’ ) ^+d^] ^ 


where dA^ = liiRdR, 0< R 5 R^ Let u= then integrating over o' gives 



1 1 


? f f (u ■♦■v a) dv du 

for a 5 (h/Ro) , with S(a) ^ ^ \ [(u+v+a)2 -4 

For the present geometry, a=l,96 and numerical evaluation gives S~0. 2710823, 

*^2 -^1 

Using Rq= 0.396875 cm, 1^ (|ig*pulse -cm -stef ) is given by a datum 
PCug/pulse) as Iq = 2.373F. 
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APPHNDIX ‘I 


si.^riin coujicroR solid angld cai£uij\tiqn 


The QCM signal results from particles scattered in a plume slice defined 
by the collimator. For each volume element in a slice, the QCM collector 
opening presents a solid angle n to the scattering source; reduction of the 
QC^l data to find the source distribution thus requires an expression for the 
source weighting factor i?/4n. An analysis of the collimator geometry follows, 
see l>igure 4-L. In the figure, the source or scattering point is It. = Cx , y , z ) 
in coordinates attached to the collimator. The collector is at an angle t relative 

to the z-tuxis (Figure 4-2) and is defined by [x] - li/2, |c| - w/2. The aperture slit 
of the collimator is defined by r=r^, 1^1 ' x=rsin0, y=rcosn. 

The collector imit nonnal vector n= Co,cosT,-sim) and a collector area element 
dxd^’= (JA is located at it^=Cx,csinT , f,cosT). The corresponding solid angle 
differential is 

1 % n • K 
dii *- — » — clA 

ir 

where it = Rc,-lt^,. This is to be integrated over the intersection or overlap of the 
opening ami the projection or mapping of the aperture slit from the source point 
to the plane of the collector. (The penumbra corresponds to a partial overlap, 
tlie umbra to a null overlap, and the illuminato to full overlap.) The mapping 
is givcni by 

I 

X =“p- [r^cosT (y^sino -x^cos(i ) +sinx (x^ z-r^ z^sinO )] 
r, coso) 

where 0 = cos t (y,.-r^coso ) -sint (zg-z). 
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Figure 4-1. Collimator Perspective 
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'rhe Jacobian is 


1^ 

^ fy COST - ZgSinT )^6cgSine + /gCose - r^.) 

3(o,z) 1) 

The mapping is nonsingular as long as D>0*. only the region y_>r is (and need be) 


considered. Since R_<< R_, the solid angle is approximately 

w 9 


ilii. 

tU 


(y^cpsi - ZgSin-r ) 

r7 


1 +j^ [^XgX + CygSiiiT+ ZgCosx) cj 


1 2 ^ 2 

with an error less than , i.e. , less than 1% for R_*r . With 

3 Rg*^ ^ ^ 

Rg» r^ for most volume elements, this expression gives negligible error. 

Let BH^c/y. <1 and x_= psincji, y = pcoscji. The inverse mapping through order 

S 5 S 

'' X \^ , is 0 =({: + (l-6cos(|)) (l+3coS(j)5)n 


Z=pcos,|>Sg+BCl-Bcos,j>) j^(l+Bcos(J)C )(| — ygCOST-^ZgCos<|') -SsCl-scostf)) n^ZgCOS(f)l 


where sin<})+ sintcos(|),n= ~ coS(ji-|- sintsincjt. 

^C c c 

'Ibe pcnumlu'u region with respect to G at 0=0^^ is for 

2 

-(I-BCOSO,,) coso^^ - ^ Sint sino + A(f ) 


X 


+ B f f Sint -C (^sim)' 


2i 


where 

A = Cl-Bcoso„ )sinG„ cosG„ 

fl ci cl 

2 7 2 

B = cos G„ -sin“0„ + 23sin 0^ cosG^ 

Ed d d 

2 

C = sinG^(coso^+ 3 sin 0 „) 

cl d d 
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Since e„ is relatively large and the region sniall, any variation in the presumed 
small source function is neglected and half the region is included in the 
illuminato, the other half in the umbra: 
dn * 0 if |<jil > 

Ixj - h/2 if |«})| - (j»* 


where (l-0cos0„) 

ti ci 


(i Sint) sin0g-A( +C ( ^ sinr)' 
c c c 


The problem of calculating n is now reduced to that of the peniunbra region with 


respect to z (or i;). Define c'= ^(z=- — ) and max(s #"w/2) , ^2= niinCcT 7 )• 

2 

Then ifO if ^2 combined penumbra- illuminate is for < x, < ?2' 

h/2 


+ w 


a * 


(y^cosT'ZgSim 



r:: 


1 


dx (c2“Cj) 
h/2 


3x 3 

1 + — g-x + ^ (ygSina- + ZgCost) (52+^1) 

R^ 2R 

s s 


The integration is to be done numerically. Angle is to be calculated with 
C replaced by equation for x, can be put into an approximate 

quadratic form: 

K2 - Ki C ♦ = 0 

K KqK2 

giving (1+ *^o “ Cz-3cos<|» ZgP)yg 

^ 1 

wliere 


p(x) l-(l-3cosi(.)r 


r 2 -I “ 

siii({i+hcoS(ii [cos(i)+g (2"3cos ({i)J r 


2 

Ki = ygCosT(l-3cosij))+sinT [z-ZgCOS())CcoS(j)+3sin (}»)] + Cl-3coS(f>)sin(f)Cosij) 

2 X 

[r cost+Ccos<})+3C1-3cos (j>)) z sim ] ^ 


2 r ^ T 

Y ^2 - Cl*3cosiji) sin (|>sinT L cosr+(l-33cos4i) cosij) ^ sinr J 

2r^ 
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i> 




* 



Each source point has as its weighting factor the collector solid 
angle sj/ 4 tt times the appropriate volume element size. 


A shnilar formulation holds for a surface source function. In particular, 
the solid angle factor may be used when weighting wall backscatter if attenuation 
is neglected. Also needed in this particular situation is the additional far '■or 
of the cosine of the angle between the QCM solid angle direction and the wai^ 
nomal. (The wall boimdary is considered to be a smooth surface defined by 
the inner exti'emities of the anechoic chamber. I The wall is elliptic with 
its center located at 4.445 cm downstreasn of the PPT exit plane and with 
semimajor axis a =^150. an and semiminor axis b = 100 cm. The wall is given 


as r. = b 
w 




The wall normal vector is pi’oportional to 


dr, 


(r^^coso, r^sino, -r^^ 


w. 


dz 


•)' in plume based coordinates, where -r 


dr. 


w 


'b\“ 


w 




dl 


Also, the wall area element dA = r,, -nr dzde^, where 1 is the axial wall arc 
’ w w OZ ’ 


dl 


length, ^ 


= b 


! V 2 /Z-Z \ 


*5 


The solid angle direction is sbnply (r^^coso, r^^sino,z) translated to QCM based 
coordinates (without rotation through dip angle, see Appendix 2). The cosine 
is calculated by taking the inner product of the normalized direction vectors. 

Total wall backscatter into a, QCM can be easily foimd by summing for 

small, discrclo increments Az,A0. 

Fortran Subroutine Listings 

The subroutines used in the upper boimd wall scattering correction estimates 
and source function .analysis follow. Subroutine names, usages, and argtmicnt 
list definitions arc: 
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NAME 


USAGE 


ARGlJ^lE^f^S 


( 1 ) UEFLC 


(2) WSF 

(3) MATOIX 

C4) IVFUN 


(5) OMEGA 


(6) COEFF 


Used to calculate corrections AIJHA - dip angle: ZA-axial position 


due to wall scattering 

Gives normal wall intensity 
Calculates influence matrix 
of volume source function 
Does coordinate transforma- 
tions and calls solid angle 
subroutine 


Calculates Q. 


Calculates coefficients 
used to find Cj 2 


of QCM pair A: IQCM-QCM number: DT, 
DZ-intervals A 0 , hi for integral sum 
R,Z - coordinates r ,z 

W 

AM-array name: M,N-row and 
colurai size 

X,Y,Z or R, THETA, Z- source 
coordinates; WDR- -r^ -gj; ALPHA, ZA, 
IQCM-gives QGM position, see above; N- 
number of intervals for Simpson's 
integration for fi; LGL-logical 
variable, if FALSE, then WFUN 
multiplies Q by cosine factor 
XS, YS, ZS - source coordinates 
relative to collimator: TAU- collector 
slant angle; H,W- collector hole size; 
ATHBTA, A- aperture opening angle 
and width; RC- collimator can radius; 
N- integration intervals 
QO, Ql, Q2, RiO, Rl, R2 coefficients; 

X- collector surface coordinates; 

Z- aperture edge position; YS- 
source coordinate 


(7) LIMITS 


C8) RM 


Calculates C 2 1 i X-collector coordinate: ZL- 

aperture edge position; YS, ZS - 
source coordinates; W2- collector 
half -width; ERR- maximum quadratic 
term in calculations of c*; ZP, ZN' 


(C2-C i). Cc2 nj) 

Calculates plione boundary THETA, Z- plume coordinates 

radius (TPA-Tan of exiaansion 

angle) 
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1 - 

?- 

M * 

‘ 4 

b* 

7 * 

) 1 
lo- 
ll^ 

1? * 
1 

1 H * 
1S< 
If * 
] 7 -» 
1 s • 

1 c i 
:’ 0 « 

?i * 

y? * 


rs 

P.L 

P7 

yo 

?■; 
-i [, 

I '1 

/.3 

7 <* 

<«. 


riMCI ION FLC <ALFHA*?A»IUCK*t>J •P^> 
DATA TINT,7If'T/Tfirt.n tir»4.<»4f /P/ino.o/ 
x=n ,0 

V=0.0 
7~n 
l.'D»? -[' 

A = W MINCXtYtZ »ALPl'A*ZA,I0CM*6,.f ALSF .) 
SUM:i().0 

IT = TINT ADT + 0.5 
IFUT.lT.l) TT=1 
A = FtOAT UT) 

IS = T INF /A 
I? =7INT/l)7 + n «,»■' 

FfCIZ.LT.l) TZ=1 
A-Flf'ATf 17) 

7S=7INT/A 
7 = 0 .‘i*7S 

i ti KZ = 1»I7 

vDRpI'.AAAAA* 7-1 .'•7556 
K= 1 .0-1 .c'iF-0«WDh*WDR 
IF C f. .LE .0.0) GO 10 pn 
Sf =I-*S(>|KTIR) 

F = 1 

IF (K.Ll.O.C) GO 10 ?() 

P = H ♦SDKT (Fo 1 
rF = SF*W';F IR» 7) 

IF'ETArr 

DO Hi KTil,IT 

A=WI 1 <ftfTHFTA,7t<DR) 

IF ( y'.LF.f’.O) GO 10 10 
SUN=SU»****SF 
10 P'E 7/' = TDfcT A+ IS 

i'l ' 7 = 7 *?S 

RTFUfrl .7AE.3P05F -? * T G* 7S * SU »>' 
t-TTIir' 
f K- D 


1 « 
?* 
r * 

A A 

•J * 


FUN n ION WSF <R .7) 
I'F.IrATANr ( F'*?-<f.A4D) 
.v$F=* .'’?F-o*rXF>(-7.7 
icflllFN 
r i\ n 


ORIGINAL PAGl K 

OF POOR QUAUTf 


- 140 - 



1 * 

£• 

3* 

4* 

f.* 

7* 

<' * 

10 * 

11 * 

1 ?* 

]l* 

34* 

i:* 

16* 

17* 

16* 

IS* 

?c* 

21* 


SUBROUTIM' 

PARAMETER NDP= 40 ♦ NSP = 3». 

REAL AM<\'DPfNSP) *2V <NDP ) *ALPHV CWDP) *R1 (NSP» *R2CNSP) « 

C TKAfSP) fT2(NSP)f21«NSP)fZ2(NSP) 

INTEGER Mt*Ur0V(NDP)*lR*lTtI2 

1 FORMAT <2I3flP3E12. 5) 

2 FOR'^AT ( 1P2E1 2.5,1 3) 

3 F0RMAT(1PFEJ2.5) 

4 PORMATaMi*4X21HR *THETA,Z 1NCREMENTS=1P3E12 .5 > 

5 FORMAT (3H(.,5X35H0P HA ALPHA IQCM/ 

C #5X13, lXlP2E12.5,‘jX12>) 

6 FORMAT (1H0,5X2HSP ,5X4HRMIN,8X4HRMAX,6X8HTHETAMIN,4XeHTHETAM AX* 

C 6X4HHM IP *fiX 4HHH AX/ C 5X13, IX 1P6E 12.5)1 

RE AD (5*1) M,N*DR#DT*D4 

READ (5,2) (HVf I)*ALPHV( 1 r- * 1 Q V C I ) * I = 1 *M) 

RE AO (5, 3) (Rl( n*R2(T ),T1 (7) •T2(I)*21(I )*22 (1 )*1=1,N) 
WRITL(6*4) DR*DT,D2 

URITC(6,5) ( I.2V(I),ALPHV(I),i9V(l)*I = l*M) 

WRITE (0*f.) ( I,R1( I),R?(1)*T1(I}*T2(I),21(I)*Z2(I) *I=1*N) 

PC iro 1=1, M 
y = T . t- 


22 * 

23* 

24* 


Y = f/.C 
7 = 5' . 
A = 


WFUP( X,Y,2,AIPHV(I),2V( I), I(3V(I),6,.TRUE.) 


25 * 

no 111 J=1,P 

26* 

suM=o. r 

27* 

A=T2(J)-T 1(J) 

2F* 

IT = A/DT*': .5 

29* 

irUT.LT.l) 1T=1 

30* 

H=rLUAT(] T) 

31* 

TS=A/D 

32* 

THE TA=T1 ( J)+C.5* 

73* 

A = H2 ( J )-H 1 (d ) 

34* 

lZ = li/L'2*f .5 

36 * 

IF(IH.LT.i) 12=1 

36* 

6=FL0AT( 12) 

37* 

7$=A/h 

76'* 

no 5i, KT=1,IT 

’S * 

2=21 <d )♦■. .5*25 

4 ( * 

CO * K2=1,I2 

41* 

A=RM ( THFT A,2 ) 

42* 

A = AMIN'l ( A ,R? (d) ) 

43* 

A=A-R 1 #d> 

44 * 

IF (A.LF.I'.C) GO ' 

4h * 

Tf>,=A/nR + '^.C 

46* 

ir(IR.LT.l) 1R=1 

47* 

0=FL0AT(Ik> 

48 ♦ 

RS=A/B 


TO 32 
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(vD* •.i'*PS 

b( • 


!'0 3. KK=i,lR 

Ul* 


1 <K »1 HLl A*?»5G,G ) 

b?* 


A=A*R*PS 



SUK=SUM*A 

54* 

31 

P=R*RS 

55* 

3? 

co^TI^ur 

5f * 

4' 

2=7*?? 

17* 

> 

THFTA=THrTA*TS 

51* * 


A=1.74532V3l -2*T?*7S 

59* 


aM<I,J)=A*SUM 

* 

If' ( 

rUNTlNUE 

Gi * 


RFIURN 

f-Z* 


F HV 


1* 

1 UNCTION WFUN( XfV ,Z»ALPHA»ZA 

2* 

RCAL X*Y,7*ALPUA*2A»?C 

3* 

INTlfTFR 1QCK»N»J 

4 • 

LOGICAL LPL 

5* 

riMfNSiPN SC(7),V(3) 

* 

C ONKON/SAC/SC 

7* 

f'ATA DAfnK,n?*RC« A*H*W,HTD/i; 

fi* 

C 1 .G35,t .5 35*1 .74532935- 

« 

X=X 

1C* 

Y = Y 

n* 

7=7 

ir.* 

CrALPHA-CA 

13* 

PL=ORF f.A<RTn*nA»RTC*C*Hf W,5v 

3 4* 

J= IMCM*1 

It* 

J = J/? 

H * 

*1=0-1 

17* 

fL=FLC AT ( J) 

Ifi* 

YT = 38. t*lF.f'*RL 

19* 

7T=i A*P7*RL 

PC * 

7T=-ZT 

71* 

*.1 = 1 fi C N / ? 

? L * 

o=incp-j-o 

r 7 ^ 

C=RTD* ALPHA 

?4* 

?=?] N( r } 

75* 

c = c'o? (c ) 

25* 

Cl' TO 1 1 

?7* 

fNTf'Y Wn (R,THETA,2,yrR) 

?{» * 

LPR = WI>R 

P9* 

PL=RTr*THf TA 

7- * 

>= R*C0?|RL) 

71 * 

Y= R *? INIRL ) 


ORIQINAL PAGE |§ 
OP POOR QUALfTY 
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32* 


( NTSV WF2CX, V,Z> 

33* 

1 f 

>S = X 

34 * 


Y5=Y*YT 

35* 


7S=7*7T 

36* 


RL = PX 

3A* 


1»^CJ.F0.1) RL = -DX 

36* 


XS=XS*RL 

39* 


AF =1 .1' 

40 • 


ir<LGL) (in TO 2C 

* 


V n )=xs 

42* 


V(2>=YS 

43 * 


V (3) =2S 

4A * 


CAt.L FL2NRMf 3* V*RL ) 

45* 


V ( l)=x 

46 * 


V( 2)=Y 

i*T* 


V(3)=UDR 

4n * 


CALL 6L2NRM< 3. V,AF 7/ 

49* 


AF = (X*XS*Y*Y6 + UPR *2S) / IRL *AF ) 

51' * 

2 ; 

CONTlNUr 

51 * 


RL=YS 

52* 


YS=VS*C*7S*S 

' ’>* 


?fr7R*C-RL*S 

14* 


R1 XOMFOA2 <XS .YS,; S ) 

55* 


WFUN= 7 .95 774 72 E -2 *RL *AF 

56* 


RE TURN 

57* 




1 * 

n * 

>=^t* 

f-. * 
t * 

u*- 

?• 

!'■ * 

n * 

13* 
14 * 
1 *: * 
If* 


rUNCT ION 0^'^;^^<XS*VS*2S♦TAU,H^U* ATHETA*A»RC»N) 

tn M XStYS,7StTAU,H*U,ATHETA»A*RC,S ,XtA2f U2 fFRP»7D,ZA 

INTPOrR N,J 

ror^MON/SAC/KCPfSINTAU tCOSTAU*BrTA»SIKPHl*COSPHI*CF 
1 FUR'^AT (4XbHBETA = lPElC,3f2X3HRS = lP3ElC.3> 

r- f OR NAT (AXRHFRR =1PF 1 0 . 3»2X3HRS = 3P3E1 0 .3) 

3 FORNAT(4X6HOMEOA=1PU C.3,2X3FRS=1P3E10*3) 

TAUR = 1.7453293E -2*TAU 
the TAR =1.745 32 931' -2*ATHFTA 
SINTAL'=SIN(TAUR) 
rOSTAU=COS(TAUR) 
fINTHT=SIN(THFTAP ) 

COFTHT =C.OR(T Hf TAR ) 

A2 = i .5*A 
h? = '-. 5*H 
V2=;.5*W 
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i 

i 


17* 

IH* 

X9* 

ac* 

22* 
23 # 
2<»# 
25# 
2fc# 
27# 
28# 
29* 
3C# 
31# 
32* 
33* 
34 * 
35# 
3fc# 
37* 
3B* 
39* 
4 (. * 
41* 
42# 
43* 
44 • 
45* 
46* 
47* 
4H* 
4'=)* 

C ^ 

*..» k " 

51* 

52* 

54* 

55* 
56* 
57* 
58* 
59 * 
f * 
61 * 
62 * 


64* 



RCP-KC 

r.!=N-N /2 

S=MOAT(N) 

fNTRY 0MCnA2<XS,YS»ZS> 

f/MCGAsf .( 

RS =XS*XK#YS *YS#ZS*ZS 
rs3sRS*S«KT<RS ) 

PHI=ATAN2 (X£,YS) 

S TNPHl =S1N(PHI ) 

COSFHI=COS(PHI ) 

IF <YS.NE,'.,in BtTA=RC/YS 

IFtC.O.LE .YS.AI^D. YS.LF.RC.ANC. ABS(XS>.Le.RC.AND.A8S f ZS) .U. A) 
C WRITEI6,1) BETA#XS,YS*2S 
IF (AbRIPHD.Cr.l.O) RETURN 

IF(YS.LE.PC.OR.APS<2S*BETA*COSPHI>.GT.I A2#-y2>> RETURN 
CFA=1 .r-PrTA*COSTHT 
F. = CFA*S1NTHT*CDSTHT 
L‘=SINTHT* (CPSTHT + BETA*SIN7HT*SINTHT) 

CF = l,r>Bt:TA*COSPHI 
F1=YS*CPSTAU-2S*S1NTAL' 

F2=YS*£IN7AU+2S*C0STAU 
C S1PPS0N*S RULf INTEGRATION 
X = -H2 

E p R = p . r 

CALL LIBITSiy, A2# YS»ZSfW2#ERR*ZD»2Al 

F 1 = 2 A 

f.2rZD 

r?=x*zo 

E4=ZA*?D 
DO 14 Irl ,N 
J= 1/2 

X=X#H/$ 

CALL L IWI T£( X*A2*YS,Z5*W2»ERR,2DtZA) 

WF = 2,L' 

IF<J.E«.1> WF=4.L 
JFU.CO.M) WF = l,t 
r 1=F1*WF*2A 
r2=E:2*wF* zn 
L3=F3*WC*X*ZD 
r 4 = F4-*WF*ZA*ZD 
1 . CONTI NUf 
WF = 3. ■•*S 
Fl=f:i /WF 
^ 2 = F2/U'F 
r 3 = C 3 /UF 
I 4 = r 4 / WF 


« 
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c 


ir<rHR ,!.T,< . n UKlTf:<f«2> ERR*XS*YSt2S 


f;.‘‘ * 

r6* 

‘ t.i* 

r.9* 

7 \' * 

71* . 

17 * 

73* 

7<t* 

75* 

7C.* 

11 * 


X-H2/RC 

WF = l‘.'?*ri *SI NTAW/HC 

WF = THF.l AR*CF A* fSTNTHT*WF*C.*kF*WF-B*X*X) 

lFfABS(PHI).(vT.UF) RTTURN 

WF=3.''/PF 

X = E: 2 *WF*Xf;*C 3 +C. 5 * WF*F 2 *EA 
WF=F 1 *H*X/RS 3 

IFFUF.LT. ' . J) V,'R1TF.(6*3) UF,XSfYSf2S 
HMEftAs AF'A ''I { WF,C. r ) 

r- r TURV 

rwn 


1 * 

c.' * 

3* 

« * 
f * 

(■ * 

r 7 * 

h * 
1 * 

,1 : * 
1 1 * 
12 * 
13 * 


MiFROUTIN^ rOf FKQCl*RU 02 *Rt,RlfR 2 »Xt 2 tYr>) 

R F AL (V 1 f C 2 » R J » M 1 » R 2 « X, 2 , V S 

rrFMorj/SA(:/Rr.sTAUfCTAU,pt ta*sphi*cphi»cf 
r:'=ys*z 

:n = CF*CTAIJ* YS* (1. C + BeTA*SPHI*CPHI*X/RC)+Z*?TAU 
0?.=CF *STAU*CT AU*SPHI*SPHI 

R = RC-CF*X* CSPH1^E.5*CPHI* (CPHI + BETA*( 2.0-3. C)*CPHI*CPHI >)*X/RC) 
F’= CPHl*Rf 

r-;l = CPHI*FETA*SPMI *SPHI-CF*SFFI*(CPHI+BE TA-3.C*BETA*CPHI*CPH I)*X/RC 
Rl=CPHI*STAU*Rl 

R2 = CF* (u.5-1 .‘j*BF TA*CPHI )*SPFI*SPHI*CPHI*STAU*STAU/RC 

PFTURN 

: "'n 


0/v 

Oa 


1 * 

FUNCTION RF<rMFTA*Z) 

2* 

FFAL THFTA.Z 

3* 

haTA KF,TFA/7.0» . .267 

h* 

t =.rf.TPA*7 

5* 

p = 6#6 0666*‘>7r —3 *Z — 2 .96 

h* 

n=l,t-F*D 

1 * 

' - IC'O.l *S'3FT tP ) 

b * 

= A*P ) 

n* 

PflURN 

1 u * 

r ND 
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1 * 

2* 

5* 


* 


7* 

fl* 

1. 

lU 

124 


134 




^ *=■ 


16 ♦ 
17* 
IH* 
ly* 
2;* 
21* 
2 2 * 
?3* 
24* 
?5* 

?t) * 

27 * 
rf>* 

pf?* 



CUt-HOUTlNf:, I IK ITS(X*ZL,YS,ZSfW2*[:RR*2P*ZN> 

RPAL X,ZLf YS*ZS*U2.ERP*ZP *ZN*Sf 00*01 #Q2»Rt*Rl*RP 
['IHENSION S(7) 

COMKON/SAC/S 

CALL COfFF (0 0*01, (?2*Rr.,Rl,R2*X*ZL*YS) 

fasOO-RL*7f5 

flsOl-Rl*ZS 

(;2=02*R2*ZS 

ZP=rO0*G2/(G: *01) 

IKAf S<zr 1 .OT.tRR ) ERRSABSCZP ) 

ZPr-Oii* <1. r^ZPl/Ql 

CALL C0tFF<0i *01*Q2*R('*Rl*R2#X*-ZL*YS) 

Cv=Q?-pn* z.s 

(.1=01 -R1*2S 
02=Q2*R2*2S 
2W=on*02/((31 *01) 

1MABS(2N).0T.ERR> ERR = ABS(ZM 
ZNsOf * (1, r*2N) /Oi 
n? = ZP*2K' 

■*=«2*lif2 

!r<ABS(Q2 ) .GT. W) 0 2 = S IGNf W *02 ) 

Z2=A^INK?P*W2> 

Z\' = APAX1 (ZN* -W2> 

Zf-'=ZP-2N 

ZP = A^"AXi < ZP* *• U) 

ir (/p,(;t*c*r>) or=zw*ZN 

ZN=ZP*Q2 

PE TURN 

FM) 


fWOlNAl PAoe 
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AI’PUNDIX 5 


EFFECTIVE QCM SLIT ANGLE 


It is desired to reduce the QO^i data (in pg/pulse) cf plume backscatter 

to a iocat intensity value. To do this, the effective collecting area-solid 

angle product needs to be calculated. For a given element of area on the QCM 

surface, tlie solid angle is defined by the aperture slit area. Coordinates 

relative to the collector have the x-axis as the tilt axis and y-axis as 

the look direction (see Figure 4-2, Appendix 4). The aperture slit is r=r , 

ji)| - -> = rdV’, | 2 | - a/2 for x » rsino, y » rcosQ. Since the QCM surface 

normal is at an angle t to the look direction, (see Figure 4-1, Appendix 4), 

on tiie collection surface y *r,siiiT, c =ccost and the area element dA = dxd^. 

The effective collection area is cosG cost dA^. For given x, let ft = \pej.tuj.e 

-11., n then (h: = R (I^'n) r, dixls where = r, “+x'^+r/+z'^ -2zccosx-2i' 

collector c c c 

(x.sinn+r,simcosu) and (R'n) = R“-v“, v“- (xcoso-tsinTsino) +(z-?cost) . 

Integrating over z from -a/2 to a/2: 

G.. a/ 2 


Jcomxhi - 


;i 

I cosi> do 

U-tnn-'f 

k 

U 


-r, 7- 


ua 

T 


I 


wdz 


wlK're 11 “ - r^,“+x“-':.“sin“! - 2r^, (xslno+f.sinrcoso) and w = 1- Cl-vVp2)'^. The 
argiuiieiU of tlio inverse tangent is small; tiic fiuict ion is ro])laced by the first 

A I 

two terms ol' its Taylor series. IVitliin an error of order — ) - 10 ‘ : 

‘ c 


w 


and 


^osodii= 


costxlo 


in: 


u +c cos x-a /^ 2r 


,„r 


(Second U'liK luis I^ec'ii reduced to lowest order.) The first tenn can he expanded 
for large r^ and r placed by 


1+ fxsinn+r,simcoso)+ ^2(^sino+(;sintcoso) 

*c 


1 ,2^ 2 J -2 

{* +'■. - T ) I 

c 


h h 


Integrating over x from ” j to ^ and over 5 from ■•w/2 to w/2, 
ycosodiklA^. = 2 I 1 W (l+c) ~ sin the small factor 


~ ~ i (rM ^a J 


;3 . 


2 r 

t) [ 


C2+cos"'0„) sin^T-1 

ci 




Since C is of order 10 , it may be neglected. The effective slit angle is 

simply =0.1 radian = 5.73°. 
c 


Intensity variations with respect to 0 (plume transverse direction) will 

be neglected asid a mean intensity as a function of the dip angle only is 

calculated. This intensity is the datiun divided by coSi cm“-rad, 

*c 

-9 7 

A^ = hw. Numerically, the divisor is 4.03225 x 1.0" enr-rad cost or 
2 

2.31 an -deg cost. 
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